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All known nuclei are made of the two nucleons, the proton and the neutron. 
Besides those two lightest baryons there exist still a couple of other stable 
(but weakly decaying) baryons, the hyperons. These were found for the first 
time in cosmic ray experiments and were entitled as 'strange' because of their 
unusual long lifetimes of about 10~ 10 sec. For their notation in the hadron zoo 
of elementary particles one introduced the new quantum number 'strangeness' 
(or 'hypercharge'). In the quark picture baryons are interpreted as being made 
out of three constituents, the quarks. Correspondingly the proton (uud) and 
neutron (udd) are made out of the two lightest quarks, the up and down 
quark. For the description of the hyperons a third flavour, the strange quark, 
was demanded. The lightest hyperon, the A-particle (uds), and the three S- 
particles (uus,uds,dds) contain one strange quark, the two more heavy El- 
particles (uss,dss) contains two strange quark and the f2-particle (sss) is solely 
^ \ made out of three strange quarks. 



Up to now strangeness remains an experimentally as theoretically rather largely 
unexplored degree of freedom in strongly interacting baryonic matter. This 
lack of investigation reflects the experimental task in producing nuclei con- 
taining (weakly decaying) strange baryons, which is conventionally limited by 
replacing one neutron (or at maximum two) by a strange A-particle in scat- 
tering experiments with pions or kaons. There exist a broad knowledge about 
single hypernuclei, i.e. nuclei, where one nucleon is substituted by a A (or E) 
by means of the exchange reaction 7r + + n — > A + K + (the K + (su) denotes 
the positively charged kaon, the lightest strange meson). Over the last two 
decades a rich phenomenology has resulted for such hypernuclei. 

However, there exist more or less no experimental insight how more than one 
hyperon behave inside a nuclei (or nuclear matter). The technical problem is 
to create within a tiny moment, smaller than the decay time of a hyperon, 
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enough hyperons and then to bring them together with nucleons to form any 
potential multi hypernuclei. By employing a relativistic shell model calcula- 
tion, which gives a rather excellent description of normal nuclei and single 
A-hypernuclei, it was found that such configurations might exist as (small or 
large) bound multi hypernuclei (MEMO - metastable exotic multihypernu- 
clear object). The reasoning for such exotic bound nuclei and the resulting 
extension of the periodic table to the strangeness degree of freedom will be 
discussed in detail in section 3. 

Hypermatter could also be realized in a completely different picture. Indeed, 
this second and much more speculative possibility was raised by physicists 
much earlier. The fundamental theory of strong interactions, quantum chromo- 
dynamics, does not forbid the principle existence of 'larger' hadronic particles, 
so called multiquark states. Today only the mesons and baryons are known in 
nature. However, there could exist states with more than three quarks. Going 
further with this speculation one comes to the conclusion that only multiquark 
states with nearly the same number of up, down and strange quarks might 
exist as stable configurations. According to very schematic model calculations 
such states could exist as metastable or even absolutely stable objects, being 
then more bound than 56 Fe. Such a very speculative form of hypermatter is 
called strange quark matter. We will start our discussion with an introduction 
to the physics of strange quark matter in section 2. 



Fig. 1. Phase diagram of (hot) nuclear matter including the strangeness degree of 
freedom - MEMOs and possibly also strangelets establish as stable multistrange 
configurations. 
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formation of either multi-hypernuclear (strange hadronic matter) objects, con- 
sisting of nucleons, A's and S's, or strangelets (strange multiquark droplets). 
To be more specific, on the average the occurrence of 20 A's, 10 S's and 1 
S's per event for Au(11.7 AGeV)Au and of 60 A's, 40 S's and 5 S's per event 
for Pb(lQ0 AGeV)Pb are expected to be created. In Fig. 1 we depict quali- 
tatively what we want to elaborate in section 4: Customarily the equation of 
state of hot and dense hadronic matter (being created in a relativistic heavy 
ion collision) is characterized by means of a phase diagram «-> T), where 
at some critical temperature and/or nonstrange baryon density eventually a 
phase transition to a deconfined quark gluon plasma (QGP) state does occur. 
However, the EOS to be passed through during a heavy ion collision incorpo- 
rates also a new degree of freedom, the net strangeness (counting a surplus 
of strange over antistrange quarks). Like the occurrence of bound nonstrange 
nuclear matter, multihypernuclear matter, or small droplets (MEMOs) of this 
new state, may be revealed. In addition, also the phase transition to the de- 
confined state is affected by the possible conglomeration of the strangeness 
degree of freedom. In particular, if the strangelet does exist in principle, it 
has to be regarded as a cold, stable and bound manifestation of that phase 
being a remnant or 'ash' of the originally hot QGP-state. We will close sec- 
tion 4 in critically emphasizing the detection possibilities of such exotic states 
by their properties and lifetimes, also in respect to the present experimental 
undertaking at Brookhaven and at CERN. 

In section 5 we finally want to give some insight on how the physics of strange 
matter can affect the physical picture of dense neutron stars. 
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2 STRANGE QUARK MATTER 



The proposal that hypothetical strange quark matter droplets ('strangelets') 
at zero temperature and in /3-equilibrium might be absolutely stable has stim- 
ulated substantial activity. Such a scenario would be of fundamental impor- 
tance, as, for example, it could explain a large fraction of the non-observable 
mass of the universe in certain cosmological models, and could modify sub- 
stantially our understanding of the structure and stability of neutron stars 
[1,2] (see also section 5). 

Some years ago we proposed that such strange quark matter droplets might 
as well play an important role in ultrarelativistic heavy-ion collisions [3-5]. 
One important goal in ultrarelativistic heavy-ion physics is the observation 
of a temporarily created quark gluon plasma. An enhanced production of 
strange particles from the QGP, especially an enhanced i^ + /7r + -ratio [6], and 
the //-^-suppression [7] were proposed as possible signatures for such a novel 
state of deconfined, strongly interacting matter. It was subsequently shown, 
however, that such indications, as indeed were experimentally found, can be 
still understood as being due to the formation of a very (energy-) dense and hot 
region of matter consisting of confined hadrons or precursor effects of QGP 
formation. Therefore, it seems that perhaps the only unambiguous way to 
detect the transient existence of a QGP might be the experimental observation 
of exotic remnants ('ashes' of the QGP), like the formation of strange quark 
matter droplets. To this interesting aspect we will turn back in section 4. 

In this first section we now summarize the reasons for the speculation of (cold) 
stable or metastable strange quark matter in the bulk. In addition, being more 
relevant for heavy ion physics, also a brief description of finite pieces of strange 
quark matter will be given (see also section 4). 

Strange quark matter or strangelets are thought to be confined (bulk) objects 
containing a large number of delocalized quarks (u...u, d...d, s...s), so-called 
multiquark droplets. Multiquark states consisting only of u- and d-quarks 
must have a mass larger than ordinary nuclei, otherwise normal nuclei would 
certainly be unstable, which, of course, is not the case. However, the situation 
is different for droplets of SQM, which would contain approximately the same 
amount of u-, d- and s-quarks. Bodmer was the first person to consider this 
new form of matter [8]; he proposed that strange multiquark clusters, being 
much more compressed than ordinary nuclei, may exist as long-lived exotic 
isomers of nuclear matter inside neutron stars. 

Chin and Kerman [9] and independently McLerran and Bjorken [10] postulated 
two reasons why such huge hadronic states should be relatively stable: 

(1) The (weak) decay of a s-quark into a d-quark would be suppressed or 
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Fig. 2. Schematic illustration of the energy levels inside a multiquark bag with two 
or three flavours. 

forbidden because the lowest single particle states are occupied. 
(2) The strange quark mass is lower than the Fermi energy of the u- or d- 
quark in such a quark droplet; the opening of a new flavour degree of 
freedom tends to lower the Fermi energy and hence also the mass of the 
strangelet (see Fig. 2). 

According to this picture, SQM should appear as a nearly neutral and massive 
state because the number of strange quarks is nearly equal to the number 
of massless up or down quarks and so the strange quarks neutralize that 
hypothetical form of nuclear matter. 

It was then later Witten who realized and hence raised the intriguing possi- 
bility that strange quark matter might in principle also be absolutely stable 
and may also provide an explanation for cold dark matter in the universe [11]. 
This would be the case if the mass of a strangelet is smaller than the mass of 
the corresponding ordinary nucleus with the same baryon number and hence 
be the true groundstate of nuclear matter! If being stable and nearly neutral, 
it could exist at all possible sizes [12], as the small Coulomb energy is not suf- 
ficient for a break up into smaller pieces [13]. Such a speculation is illustrated 
in Fig. 3 [14]. 

Presently such a highly speculative scenario cannot be ruled out. Normal nuclei 
could only transform into these novel states by a higher order weak decay; the 
decay into one single A-particle would be energetically forbidden. The time- 
scale for such collective decays is extremely large; accordingly, normal nuclei 
would remain intact. 

On the other hand, it is also conceivable that the mass per baryon of a 
strange droplet is lower than the mass of the strange A-baryon, but larger 
than the nucleon mass. The droplet is then in a metastable state, it cannot 
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Fig. 3. The possible places where to find stable strange quark matter. If being 
absolutely stable SQM might exist in forms of small nuggets (strangelets) being not 
(much) larger than normal nuclei. It might also span as a charge neutral state the 
empty 'nuclear desert' [12] within the range from Ab ~ 300 up to sizes of neutron 
stars Ab ~ 10 56 . 



decay (strongly) into A's [9]. 

For a very simple quantification of these ideas one typically models quark 
matter in bulk by a Fermi gas equation of state of interacting quarks which to 
first order in a c = g 2 /167T is given by [15] 
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Here mj and /jj denote the (current) mass and the chemical potential, re- 
spectively, of the quark flavour i=u,d,s. For the total potential the vacuum 
excitation energy BV has to be added, which corresponds to the energy dif- 
ference between the 'false', perturbative vacuum inside the 'bag' and the true 
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vacuum at the outside, and which, in return, confines the quarks, i.e. 



Q(n q ,fjL s ',m a ;a c ) = + BV . 



(2) 



m u pa m d pa and \i u = fid has been implicitly taken for an isospin symmetric 
situation. From this expression (2) the energy per baryon in the groundstate 
can be readily obtained by varying E/A with respect to the baryon density 
to describe the system at zero pressure. To fix the amount of net strangeness 
in the system the strangeness fraction f s is introduced as the ratio of net 
strangeness and total baryon number, i.e. 



Fig. 4. Energy per baryon for strange quark matter at zero temperature and zero 
pressure, as a function of the strangeness fraction f s . The dashed line defines the 
corresponding mass of a hyperonic matter ground-state. 

The energy per baryon of quark matter in the groundstate thus depends solely 
on the explicit strangeness fraction f s and on the three intrinsic parameter, 
the strange (current) mass m s , the coupling constant a c and the bag constant 
B l l 4 . For most purposes to come a non-zero a c can be 'absorbed' by a slight 
increase in the bag parameter. In Fig. 4, E/A is depicted as a function of f s 
and vanishing coupling constant, but for various choices of the bag constant 
and the strange quark mass. The minimum in the energy is achieved for a 
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strangeness fraction f s ~ 0.7, if the more popular and smaller mass m s = 150 
MeV is taken. According to the formula Z/A = (1 — / s )/2 smaller pieces of 
strange quark matter would then be slightly positive at its most stable posi- 
tion. A necessary condition for the stability of strange quark matter against 
strong decay is that its energy per baryon must be smaller than that of the cor- 
responding hyperonic matter (dashed-dotted curve), which can be determined 
as 

m H yp = f s m A + (1 - f s )m N - e B ; < f s < 1 

(fs - l)mz + (2 - f s )m A - e B ; 1 < f s < 2 (4) 
(f a - 2)m n + (3 - f s )m s - e B ; 2 < f a < 3 

where A, S and Q are the masses of the strange hyperons and e B is the binding 
energy per nucleon, which is, for simplicity, taken to be the infinite nuclear 
matter parameter of 16 MeV. 

If popular parameters within bag models, like the MIT bag model employed 
here, are extrapolated to strange quark droplets, both of the above consid- 
ered pictures turn out to be possible. However, metastable strangelets relax 
somewhat the stringent conditions on the choice of parameters of the bag 
model required by an absolutely stable state [13], namely small Bag constants 
B 3 ~ 150 MeV (which, within the same model extrapolated to finite temper- 
atures, would give a too low critical transition temperatures T c ~ 100 MeV 
for the onset of deconfinement). For large pieces of absolutely stable strange 
matter one has to include the Coulomb effect due to the still tiny net positive 
charge [13]. These have been neglected in the present figure 4. It turns out that 
then the minimum of the energy per baryon number will be shifted smoothly 
to f s — > 1 instead of f s ~ 0.7, so that the strange matter piece still carries a 
very small positive charge. 

On the other hand, from Fig. 4 it follows that for bag parameters B 1 / 4 lower 
than 190 MeV strange quark droplets can only decay via weak interactions, 
i.e., they would be metastable. The baryon density of a strangelet would be 
about 1.5 — 2 of normal nuclear matter density. For larger B-values (large) 
strangelets are instable. 

The idea of the existence of metastable strangelets was actually raised ear- 
lier than the possibility of a new groundstate of normal nuclear matter. Chin 
and Kerman [9] had employed the original favoured MIT bag parameters of 
B 1 ^ = 145 MeV, g*/4ir = 2.2 (a c = 0.55) and a strange quark mass of 
about m s = 279 MeV [16,17]. The minimum they find is shifted to a higher 
strangeness content. In quark matter one-gluon exchange is repulsive, if the 
quarks are massless and relativistic, and attractive if the quarks are mas- 
sive and behave already nonrelativistic. Accordingly, metastable strange quark 
matter in the ground state might as well be slightly negatively charged. (This 
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cannot be the case if it would exist in absolutely stable form, for than, normal, 
positively charged nuclei would be attracted and absorbed at the surface. Such 
catastrophic consequences have to be excluded.) One should remark, however, 
that the linear expansion to obtain the expression for E/A used in Ref. [9] is 
not valid, as the energy density as well as the baryon density turn negative 
for a c > 7r/8, so that a perturbative expansion is meaningless and signals 
the breakdown of perturbation theory. This basically sheds some light how 
'seriously' the calculations, the conclusions and the speculations have to be 
taken! 

Originally (strange) quark matter in bulk was thought to exist only in the 
interior of neutron stars where the pressure is as high that the neutron matter 
melts into its quark substructure [18,19,15] (see also section 5). In Witten's 
prospective scenario hot strange quark matter nuggets could have condensed 
out of the deconfined phase during the phase transition in the expanding and 
cooling early universe [11]. These would carry most of the tiny surplus in 
baryon number of the whole universe. The baryon number would remain in- 
side if the heat and entropy of the nuggets is carried away mainly by neutrino 
emission instead of mesons and especially baryons. If an absolutely stable 
groundstate would exist, i.e. if cold strange quark matter is the true ground- 
state of nuclear matter, the hot nuggets would further cool and instead of 
suffering a complete hadronization they might settle into these new states and 
hence could resolve the dark matter problem. Since then the idea of absolute 
stability has stimulated a lot of work on potential consequences in astrophysics 
[2,20]. 

Nuclear collisions, however, do allow also the detection and study of only short 
lived, metastable strangelets, which will decay by flavour changing weak pro- 
cesses. Most of the experiments are sensitive to strangelets with lifetimes up 
to 1CT 8 — 10~ 9 sec. In addition, the strangelets to be expected will be small. 
Hence, in the last part of this subsection, we already address the potential im- 
portance of finite size effects. Later in section 4 we come back to this analysis 
in somewhat more detail in order to try to pin down possible short-lived or 
long-lived candidates for the present and future experimental undertakings. 
Due to these effects and the wider range of employed phenomenological pa- 
rameters, smaller strangelets, if they exist, are more likely to be metastable 
than being absolutely stable. 

Small pieces of (strange) quark matter are described within a spherical MIT- 
bag by filling up the bag with exact single-particle Dirac-states [13,21,4]. These 
calculations have been done with a c = 0. The first MIT boundary condition, 
which ensures confinement, is basically a Bogolyubov-type condition for the 
scalar mass of the quark, i.e. M(r > Rb ag ) — > oo, and leads to the following 
equation for the single particle energies: 
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Here mi is the quark mass (i=u,d,s), with m u = nid = MeV and m s = 150 
MeV. R defines the radius of the bag, k is the angular momentum quantum 
number and a labels the eigenvalues in this quantum state k. The total energy 
of the finite system is obtained by summing the 'lowest' energy eigenvalues and 
adding the phenomenological bag energy BV for the perturbative vacuum: 

E total = E E N K Ei >a + b(-R 3 ) , (7) 

i=u,d,s K;a \ O / 



where the degeneracy N K = Q\k\. In addition, a center of momentum correction 
and a Coulomb term ((N — 2N S ) / R) 2 /60a is included. (In the original MIT 
model [16] a zero-point energy term of the form —Z /R is also included as a 
more or less phenomenological parameter for fitting the light hadron spectra. 
Such an additional attractive term may have some noticeable effect on very 
small strangelets.) The volume V or the radius R, respectively, has to be chosen 
in a way that the total pressure of the quarks inside the bag compensates the 
vacuum pressure B, or, in other words, dE tota i/dR = 0. 

Fig. 5 shows the dependence of the energy per baryon of a giant bag on the 
baryon number for = 145 MeV and for a (constant) strangeness fraction 
of / s = 1, i.e. SU(3) symmetric and charge neutral matter. One clearly sees a 
pronounced shell structure. Especially the first magic number shows up nicely 
at a baryon number of 6, that is, if all quarks completely fill the IS1/2 ground 
state (the so called 'quark- alphas' [22]). Generally speaking, the energy per 
baryon for a smaller strangelet with a baryon number N B < 50 is about 50-100 
MeV/N larger as compared to the infinite quark matter calculations. This can 
be interpreted in terms of a surface correction [13]. If one varies the number 
of up and down and strange quarks, strangelets do have a large variety of 
nearly degenerate states. This fact is interesting for the heavy ion experiments 
because these nearly degenerate states will be stable against strong decay and 
thus would all be detectable if their lifetime against weak decay is long enough 
(see section 4). 

We close this section with some cautionary remarks: 
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Bag with massive quarks 




Baryonnumber A 



Fig. 5. E/A versus A in a hadronic (MIT-)Bag model for noninteracting quarks. 
A constant strangeness fraction f 8 =l, i.e. N u = = N s , is assumed (m s = 150 
MeV, S3 = 145 MeV). 

• The existence of strange quark matter cannot rigorously be predicted by 
theory, although some general and appealing arguments for its stability do 
exist. A detailed understanding cannot be achieved by employing more or 
less relatively simple, phenomenological models. Ultimately only the exper- 
imentalists (in astrophysics or heavy ion physics) may prove its existence or 
nonexistence. 

• For smaller strangelets, as the above calculations suggest, finite size effects 
and a profound shell structure might be quite important. For very small 
strangelets the colour magnetic part of the gluon exchange should be impor- 
tant and should further favour SU(3) symmetric multiquark states to exist. 
In this respect the prominent H-Dibaryon proposed by Jaffe [23] would be 
the smallest strangelet. 

• Although the idea of absolute stable pieces of strange quark matter is in- 
triguing in its consequences for astrophysics, metastable strangelets only 
need to be smaller in mass then the corresponding hyperonic multi-baryon 
configuration, which is much higher than that of a normal nuclei. Because 
of this argument it might well be that (small) strangelets, if they do exist, 
are metastable and not absolutely stable (see section 4). 
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3 STRANGE HADRONIC MATTER 



Strange hadronic matter consists of baryons, i.e. nucleons and hyperons as 
quasiparticles are the basic constituents. In this case, there exists some knowl- 
edge about the properties of hyperons inside nuclei from hypernuclear physics 
up to two units of strangeness. They extend the chart of nuclides to a third 
dimension in addition to baryon number and isospin. The extension to it will 
open new perspectives for nuclear physics. This is what we are going to discuss 
in the following subsections. 

First we will give a short introduction to hypernuclear physics. Then we list 
all possible metastable combinations of nucleons and hyperons on the basis of 
symmetry considerations which are stable against strong interactions. If they 
are bound, they will have lifetimes of the order of the weak decay lifetime of 
the hyperons of 10~ 10 s. 

The extrapolation to systems with a large amount of strangeness will rely 
then on several basic features: the knowledge that the A-nucleon and the AA 
interaction is attractive, the existence of metastable combinations of nucleons 
and hyperons or hyperons alone, and the effect of the Pauli-blocking in the 
hyperon world. 

The properties of strange hadronic matter are then discussed on the basis of 
the relativistic mean field model, which is able to reproduce fairly well the 
properties of nuclei and hypernuclei. Here we will focus on the properties of 
superheavy elements. As the hyperons £~ and H~ are negatively charged, a 
similar situation as for strangelets occur and superheavy elements with A > 
250 and a low charge are possible. 

3. 1 Hypernuclei 

Our intention is not to give a detailed overview about hypernuclear physics, 
as there exist excellent reviews in the literature [24,25], but to give a short 
summary of the main features. 

Hypernuclei are known since 1953 as Danysz and Pniewski observe the first 
A-hypernuclei in a cosmic ray emulsion experiment [26] (in some of the older 
publications hypernuclei are called hyperfragments). In the following decades 
one has found A-hypernuclei up to a maximum mass number of A = 15 which 
is due to the resolution of emulsion experiments. The lightest one is the hy- 
pernucleus \H, a system of one proton, one neutron and one A. Another in- 
teresting fact is the existence of the hypernucleus ^Be: the presence of the 
A glues together and stabilizes the nuclear core 8 Be which normally decays 
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Fig. 6. The chart of A hypernuclei as well as S, E and double A hypernuclei known 
so far (taken from [25]). 

to two a-particles! Other examples for the enhanced stability of hypernuclear 
systems compared to their nuclear cores are the hypernuclei ^Li and ^B. We 
will present even more impressive ones in the forthcoming sections. 

In the 70's programs started at the AGS [27] and at CERN [28] for detecting 
A-hypernuclei by using meson beams and spectroscopy techniques. Here one 
uses the reaction: 



K~ + n 



A + TC' 



A special feature of this reaction is, that there exist a magic momentum of 
the ingoing kaon where the A is produced at rest inside the nucleus facili- 
tating the formation of a bound hypernuclear system. Another generation of 
hypernuclear spectroscopy uses the reverse reaction: 



7T 



A + K + . 



(9) 



This reaction yields a lower resolution, but deep lying states can be detected 
more easily [29]. The study of these reactions started in the 80's at the AGS 
[30,31] and at the KEK in Japan 



The hypernuclear chart with the hypernuclear states is shown in Fig. 6. Ex- 
perimentally one has seen hypernuclei up to a baryon number of A = 209 
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Fig. 7. The potentials and the single particle levels for the hypernucleus [ O in the 
RMF model. 

(Tungsten) where one neutron is replaced by a A. Moreover one is able to ex- 
tract the single particle binding energies from excited states (like ITa) to very 
deep lying states Isa which opens the possibility to study deeply bound probes 
inside the nucleus over a wide range of mass number. This was an excellent 
test (and success) of the shell model [33] . An essential feature of the A-nucleon 
interaction is that one has seen no spin-orbit splitting experimentally, which 
leads to the statement that the A behaves as a spinless neutron [28,34] (a 
theoretical discussion of this topic can be found in [35]). The potential depth 
of the A in nuclear matter can be derived from the experimental data to be 



which has to be compared with the corresponding nuclear depth of about 60 
MeV. Calculations using Hartree-Fock models [36] and Relativistic Mean-Field 
(RMF) models [37-42] are able to reproduce very nicely the observed trend of 
the single particle energy with the mass number. 

Fig. 7 shows as an example a calculation for the hypernucleus A 7 0. The poten- 
tial for the A is shallower than for the nucleon but of a similar radius. The Is 
state of the A is a separate state within this A potential. The surface thickness 
is bigger for the A density distribution than for the nucleons. This will cause 
a A-halo as the A is not as deeply bound as the nucleons. The single particle 
energies for the lp states of the A are very close to each other, resembling the 
weak spin-orbit force of the A-nucleon interaction. 
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The lifetime of the hypernuclei has been also measured. They decay weakly 
by emission of a pion (mesonic decay channel): 

A^P + tt~ and A -> n + tt° (11) 

like the free A releasing 38 and 41 MeV, respectively. Moreover a new reaction, 
the nonmesonic decay, opens in the nuclear medium 

A + p — > n + p and A + n — > n + n (12) 

which releases an energy of about 177 MeV. The mesonic channel is domi- 
nant for very light systems but it is negligible for heavier hypernuclei like the 
hypernucleus ^ 2 C [43]. This results from the dominant absorption of the pion 
inside the nucleus and of the larger phase space for the nonmesonic decay. 
Finally, the lifetime of A-hypernuclei has been determined to be at the order 
of 10~ 10 s even for A 200 [44], which is close to the lifetime of the A in free 
space. 

There exists also some scarce information about E- and S- hypernuclei. The 
experimental situation concerning the existence of bound E-hypernuclear sys- 
tems is still controversial. The strong process E + N — > A + N smears out 
possible peaks in the meson spectra. Nevertheless, one has found surprisingly 
small structures in hypernuclear spectra. All these structures are lying in the 
continuum. Recent analysis of the data about E~ -atoms indicate that the 
isoscalar potential depth of a E in a bath of nucleons in its groundstate is re- 
pulsive [45] which would prevent a formation of bound E hypernuclear states. 
The isovector potential can be attractive, as for example in light systems like 
4Hes+ and might form a bound state. Indications for such a bound state have 
been seen in [46]. 

S-hypernuclei have already been found shortly after the first observation of a 
A-hypernucleus in emulsion experiments using kaon beams. There exist seven 
emulsion events in the literature (for a review see [47]) showing that the strong 
conversion process S + N — > A + A still allows for bound systems. The po- 
tential depth of a S in nuclear matter has been derived by using Wood-Saxon 
potentials to be 

ui N) =21-24 MeV (13) 

whereas relativistic mean field calculations get a little higher value of = 
28 MeV [48]. 

There exist also three AA-hypernuclear events by caption of a H~ and a de- 
tection of their mesonic decays (for an excellent review see [49]). In particular 
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Hypernucleus 


±>AA|MeVJ 






10.9 ±0.6 


4.7 ±0.6 


i°ABe 


17.7 ±0.4 


4.3 ± 0.4 


13 R 

AA B 


27.5 ±0.7 


4.8 ±0.7 



Table 1 

The two observables AB AA and Baa of the known double A-hypernuclei. 

these are the hypernuclei ^He [50], \° A Be [51], and \\B [52,53]. These data 
constitute our only (and poor) knowledge about the hyperon-hyperon inter- 
actions. One defines the following two quantities: 

B AA (i A Z) = B A ( AA Z) + b a( a ~ 1 Z) 
AB AA (i A Z) = B A (i A Z) - B A (i- l Z) (14) 

where -Baa is the sum of the binding energies of the two A's and AB AA is a 
measure of the AA interaction. The experimental values are listed in table 1. 
The AA interaction is attractive and relatively strong. The value of 4 — 5 MeV 
has to be contrasted with the corresponding values of 6 — 7 MeV for the NN 
interaction and 2 — 3 MeV for the AN interaction. 

Hence, the AA interaction is much stronger than the AN interaction and about 
3/4 of the NN interaction. This has to be taken into account for hyperon-rich 
systems. 

3.2 Classification of Strange Hadronic Matter 

The strong interaction conserves strangeness and charge and the baryon num- 
ber. Therefore combinations of nucleons and hyperons, which are stable against 
weak interactions, so called metastable combinations, have to be classified ac- 
cording to these quantum numbers. 

First, in a nuclear environment, strong reactions in the medium among differ- 
ent baryon species can occur. There exist two strong reactions, which have a 
considerable lower mass difference (or Q- value) as the others: S~ +p — > A + A 
and S°±n — > A±A, releasing 28 and 23 MeV energy, respectively. Indeed these 
two reactions will play a major role in the discussion of bound strange hadronic 
systems and we will keep them in mind. All other reactions have a rather 
large Q-value of more than 50 MeV which can not be overcome by binding 
energy differences (remember that the maximum binding of a nucleon is about 
60 MeV). Nevertheless, the final products of all reactions constitutes the com- 
binations of two different species which are for a given charge and strangeness 
the deepest lying state and which are therefore stable against strong interac- 
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Table 2 

The possible configurations of metastable partners of nucleons and hyperons clas- 
sified according to their total strangeness S and charge number Z. 
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Table 3 

The possible configurations of metastable combinations of three different baryon 
species classified according to their total strangeness S and charge number Z. 

tions. All possible metastable combinations of two different species are listed in 
table 2. Next consider a composite system of an arbitrary number of nucleons 
and hyperons. One can easily see, that the stability against strong interactions 
does not depend on the number of baryons, e.g. on the number of neutrons or 
£~'s for the combination {n, Therefore our considerations of metastabil- 
ity are valid for an arbitrary number of respective baryons. Note that we only 
study the stability against strong interactions here. Whether the combination 
is indeed bound or not will be discussed afterwards. 

Next one can think of metastable combinations of three or even more different 
species. The former ones are listed in table 3. Of special interest is that systems 
consisting of {A, S°, H~}-hyperons are metastable (they form purely hyperonic 
matter to be discussed later). More than three different baryon species can not 
form a metastable combination in free space as one can see from this table. 
For a metastable object of four different baryon species every combination of 
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three of them must be one of the metastable triplets given in table 3. As this 
can not be the case there does not exist metastable combinations of more than 
three different baryons in free space. The situation changes if one introduces 
interactions which will be presented in the following section. 



3.3 Relativistic Mean-Field Model 

In the preceding sections we have discussed the input for a theoretical de- 
scription of strange hadronic matter. In the following we present an extended 
relativistic mean field model which is able to reproduce fairly well the ex- 
perimental data of nuclei and hypernuclei and implements all hyperons. For 
reviews about this model for finite nuclei see [54,55]. 

The basic idea of this approach is the meson-exchange picture, where baryons 
interact via the exchange of mesons, usually a scalar meson o which parametrizes 
effectively the two-pion exchange [56] , the vector meson u (hereafter the field 
V/j) and the isovector meson p (field R^) for isotopic trends. One starts from 
the Lagrangian 

L = Learyon + Lmcsoii + Lcoupling + Lcoulomb (15) 

with free terms for the baryon and meson fields 

L B aryon =^^(^7^ - m B )^ B (16) 
B 

Lmcsoh = \d»od,o - U(a) - -G^G, U + ^ V»V, (17) 

-\B^B, U + lm 2 p R% (18) 

where the sum runs over all baryons of the baryon octet (p,n,A, S + , S°, £ _ , S°, 
S~). The term U(a) stands for the scalar selfinteraction 

U(a) = \mlo> + V + V (19) 

which has been introduced by Boguta and Bodmer [57] to get a correct com- 
pressibility of nuclear matter. Another stabilized functional form has been 
given by Reinhard [58]. The interaction is introduced by a minimal coupling 
of meson fields and baryon bilinear forms 
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Lcoupling = ~J2 9oB^B^B° (20) 
B 

-£<fc,B*B7"*B^ (21) 



B 

-J29pB^b^t b ^ b R^ . (22) 

B 

For finite nuclei one takes also account of the coulomb force 

Lcouiomb = —F^F^ - J2 q B e^Bl^ B A^ (23) 

4 B 



where q B is the charge number of the baryon. We have used here the following 
notation 



W v = d^R u - d v W + iW x R" 
The equations of motions can be derived in the standard way 



and one gets Klein-Gordon and Proca equations for the meson fields with 
source terms coming from the baryon fields 



(d^ + U'(a)) a = -J29,B^B^B (26) 

B 

0"GW + m 2 u Vp = £ g U B*Blfl B (27) 

B 

PRp, + m% = £ g P B^BTBl^B (28) 



where 



U'(a) = ^-U(a) . (29) 



The equation for the electromagnetic field reads 

VF^ + mlAp = q B e*Blv$fB • (30) 
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The Dirac equation for each baryon 

[i-y^ + m B + g aB a + g uB l^ + g pB ^r B R p + q B e^A p ) ^/ B = (31) 

contains now potential terms from the meson and coulomb fields. The scalar 
field shifts the mass to the effective mass defined as 

m* B =m B + g aB a . (32) 

Analogously, the vector fields shift the energy and momentum of the baryon. 
The equations (26)-(30) and (31) constitute a set of coupled differential equa- 
tions which have to be simplified in order to solve them. In the following we will 
use the relativistic mean field (RMF) approximation, where the meson fields 
(and the Coulomb field) are replaced by their classical expectation values 

a - (a) , - (V,) , R, - (ll,) , A, - (A„) (33) 

i.e. quantum fluctuations of these fields are neglected. Furthermore the many 
body wavefunctions of the baryons are treated on the Hartree- level. They are 
approximated by a sum of single particle wave functions 

*b = £<bVVb , (34) 

a 

where the sum runs only over the occupied states. The contribution coming 
from the Dirac-sea is therefore also neglected (this is the no-sea approxima- 
tion). Note that in the so called relativistic Hartree-approximation (RHA) the 
quantum fluctuations coming from the Dirac-sea are implemented by renor- 
malizing the Lagrangian [56] which we will not discuss here. It is important 
to know that all non-diagonal coupling terms vanish in the relativistic mean 
field approximation or at the Hartree-level. They contribute at the next or- 
der which takes into account the Fock-term (this is the Dirac-Hartree-Fock 
approximation). The coupling term to the pseudoscalar fields (pion and kaon 
fields) and pseudovector fields vanishes in the RMF approximation due to the 
7 5 matrix which mixes upper and lower components of the wave functions. 
Also only the third component of the isovector field is present in the RMF 
approximation. 

For spherical and static systems (nuclei in its groundstate) the fields do only 
depend on the radius 

<7 = <7(|r1), V, = 5, V (\r\), ^ = W»^o,o(M) , K = <VVH)(35) 
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and the spatial components of the vector fields V$, Ri and A,- L are zero. Hence 
one gets the following radial equations for the meson and Coulomb fields 



[~^- r + U '^) ^ = ~ E 9*BP.Ar) (36) 
fid 2 \ 

yrdr 2 ' ' r + m "J V o( r ) = J29u>BPv,B( r ) ( 37 ) 
' 1 d 2 \ 

^rdr 2 " ' r + m lJ R o,o( r ) = J29pBTo,Bpiso,B(r) (38) 

~l^ r • A °( r )) = E ^ep,, B (r) (39) 



where the scalar, vector, and isovector densities are given by the sum over the 
occupied single particle wave functions 



N 

PsA r ) =^b^b = E u a ,Bvt B (rhoVaAr) (40) 

a=l 

N 



Pv,b(t) = ^bIo^b = E ^,b<bMv^,bM (41) 

Q=l 

V 

Piso,B(r) ="^B7oro,B , I'B = Yl UJ a , B tPa,B( r ) T 0,B<PaA r ) i 42 ) 



Q=l 



with the occupation number u^s for each baryon species and to,b stands for 
the isospin of the baryon. The wavefunctions for the baryons can be separated 
into upper and lower components by the ansatz 



( iG a Ar) Y \ 



F a A r ) O-T 



(43) 



with the spherical harmonics Yj\ m and the normalization 

oo 

J{\G a Ar)\ 2 + \F a Ar)\ 2 }dr = l . (44) 
o 

The Dirac equations can then be written in terms of these two components 



e a , B G a Ar) = (~ + ^f] F a Ar) (45) 



21 



+ (m B + goBcr(r) + g uB V (r) + g P BT ,BRo,o( r ) + <lBeA (r)) G a , B (r) 



and 



e a , B F a , B {r) = (+^ + ^jG a , B (r) 

- (m B + g a Bcr(r) - gu B V (r) - g pB T , B R 0fi (r) - q B eA (r)) F^ B (r) 



(46) 



with 



K a,B — ' 



(,3 a 2 ) ' Jet la 2 
, "I" 2) ' ja = la ~ 2 



(47) 



The densities can be rewritten to 



1 



N 



PsMr) = ^E w«,B(2ja + l)[G a ,B(r) 2 - F Q , B (r) 2 ] 
^(0 = TTi E <"«,*(2j« + l)[G a , B (r) 2 + F a , B (r) 2 ] 

1 * 

Piso,B(r) = ^B^ja + l)r 0)B [G a , B (r) 2 + F a>B (r) 2 } 

a=l 

The occupation numbers are 



+ 1 



(48) 
(49) 
(50) 

(51) 



for closed shells. Otherwise one chooses the occupation number in the interval 
< u) a ,B < 1 according to a schematic pairing 



v a ,B - ~ 



£a,B ~ ^Fermi,S 



\J (ta,B — ^Fermi,_B) 2 + A 2 



(52) 



analogous to Skyrme-Hartree-Fock calculations [59] where 

U.2MeV 

A = 



(53) 



is a constant gap derived from nuclear pairing properties. The Fermi-energy 
£Fermi,B is evaluated in such a way, that the baryon number is conserved for 
each species. 
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Fig. 8. The single particle energy of several hypernuclei in comparison with RMF 
calculations. 

The coupled equations (36)- (39) and (45), (46) together with the above ex- 
pressions for the densities can be solved iteratively using standard techniques 
(for details see [60]). The parameters of the model are fitted to the properties 
of spherical nuclei [61,60] where we will use the stabilized parameter set PL-Z 
[58] in the following. The RMF model gives a description of nuclei as good 
as the conventional Skyrme-Hartree-Fock calculations [55]. Also the binding 
energy of hypernuclei can be nicely reproduced (see [37-39], for fits to recent 
hypernuclear data see [41,42]). 

Fig. 8 shows the measured single particle energy for different hypernuclei in 
comparison with RMF calculations. One sees nicely the different single particle 
levels of the A. Note that the spin-orbit splitting for A levels is negligible 
therefore the shells can be labeled by s-shell, p-shell and so on. Bulk matter 
(A — ► oo) corresponds to the crossing of the energy levels with the y-axis. All 
levels merge to the value of the A potential depth in nuclear matter of about 
27MeV. 

It turns out that the two coupling constants of the A (g a \ and g w \) are strongly 
correlated because they are fixed by the potential depth of the A 

Ua = -g a Acr - PcjaVo (54) 



in saturated nuclear matter [62]. Hence one can choose for example SU(6)- 
symmetry for the vector coupling constants (see table 4) and fixes the scalar 
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coupling constants to the potential depth of the corresponding hyperon. 



Nevertheless, this model, hereafter named model 1, is not able to reproduce 
the observed strongly attractive AA interaction as can be seen in fig. 9 irre- 
spectively of the chosen vector coupling constant. Model 1 does not give the 
necessary strong attraction between the two A's. An additional interaction has 
to be invoked between hyperons. This can be done by using a SU(3) symmetric 
model and introducing two additional meson fields, the scalar meson /o(975) 
(denoted as a* in the following) and the vector meson 0(1020) [63,48] which 
we will discuss in the following. 



3.4 Hyperon- Hyperon interactions 

The additional Lagrangian for the new mesons is straight forward to derive 

L = L + L^ son + Lc^piing (55) 

with 



LmLh = \ (<V W - m>* 2 ) - + \m\^ (56) 

^Coupling = - E 9a*B*B*B<T* - ]T g^Bl^B^ (57) 



B 

and the notation 



Sua, = d^v - . (58) 

The new meson fields give new potential terms for the Dirac equations 

£ a ,Bf a ,B(r) = 7o • V + m B + g a Bcr{r) + g a *Bcr*(r) (59) 
+g^Bl°V (r) + g4,B^°<f>o(r) + g P Bl°Ro,o(r) + q B cy° A (r)^(p a , B (r) 
and two new meson field equations 

(-A + ml) a*(r) = - ]T g c * B p s , B {r) (60) 

B 

(-k + rnf)4>Q{r) = ^g 4>B p v , B (r) (61) 
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Table 4 ~ ~ ~~ ' 

The vector coupling constants in SU(6)-symmetry relative to gNNp- 

which are written here in the static RMF approximation. 

The vector coupling constants to the 0-field are given by SU(6)-symmetry (see 
table 4). How does one get these relations? 

Here we discuss the coupling scheme for the baryons to mesons in SU(3) 
symmetry. The baryons are members of an octet, the mesons of a nonet (a 
octet plus a singlet) under SU(3) symmetry. Coupling two octets (for example 
two baryons) results in 

8® 8 = 1©8 1 ©8 2 ©10©10*©27 . (62) 



These terms have to be multiplied with the octet of the mesons now and the 
result must be an overall singlet, so that the interaction (and the Lagrangian) 
is invariant under SU(3) symmetry. As there appears two octets in the above 
product, there are two possibilities to form an overall singlet in the Lagrangian. 
They are usually called symmetric and antisymmetric coupling. The latter one 
is given by 

L / = ~g% ■ ifimn^m^n<Pi , (63) 



with the antisymmetric structure constants of SU(3) 



fimn = -^Tr([A m , \ n ]Xi) (64) 
4z 



and therefore also called F-type coupling. The octet coupling constant g$ is 
a universal coupling strength. One can rewrite the structure of the above 
interaction in matrix form. We normalize the SU(3) matrices as 

TiXiXj = 25 i:j . (65) 



and define the baryon and meson matrices as 



B = -^=A m ^ m B = -^=\ n ^) n M = (66) 



25 



This gives the compact notation 



L f = -g 8 -Tr([B,B]M) (67) 



The case for the symmetric coupling is analogous. Here the symmetric struc- 
ture constants defined as 

dimn = ^Tr({A m , A n }Aj) (68) 



give 

L d = -g 8 ■ d im „iM>„& = -g 8 ■ Tr({B, B}M) . (69) 



This is also called the D-type coupling. 
The baryon matrix is given explicitly as 





( i yo i i A° 


£+ p ^ 


B = 




-75 S ° + 7a A ° » 




V ~ s 


WO 2 AO 

7e A / 



(70) 



and analogous for the antibaryon matrix. The general form of the Yukawa 
coupling is a mixture of D-type and F-type coupling 

L sv{3) = -g s aTT([B,B}M)+g 8 (l-a)TT({B,B}M) . (71) 



Here a denotes the F/ (F+D) ratio. For the vector mesons one gets the follow- 
ing coupling terms 



Lsu(3) = -gNNpNrNp - #£A P (£A + A£)p - g^ErEp - g SSp ETEp 
-gNNuj s NNu 8 - g AAuJs AAuj 8 - sf S Ea, s EEu; 8 - #ssu> 8 _SSu; 8 _ 
—9ank* (NtAK* +AtNK*) - gxNK^NrZK* + HtNK*) 
-Qeak* (ArEK* + EtAK*) - g~,YiK* (prEK* + SrEiT) (72) 
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and the meson-baryon coupling constants are now related to the F/(F+D) 
ratio and the overall octet coupling strength by 

9NNp = 98 <?EAp = TJS^I — «ps) fi'SEp = ^gotps 
gNNLUs = ^fl'8(4Q;p S - 1) fi'AAws — — ^^8(1 ~ OL ps ) #EEw 8 = ^^(1 ~ 

#ssu, 8 = -^s(l + 2a ps ) Qknk* = -^s(l + 2a ps ) fi^A^* = 5(1 - 2a ps ) 
^sAx* = ^^ 8 (4o!p S - 1) g~ZK* = -gs g~~p = -g(i - 2a ps ). 

As noted before, there is also a singlet state for the mesons which couples to 
all baryons with the same strength 

Li = —gi(NN + AA + EE + EE)ui , (73) 

which means that 

9nnuh = 9aau>! = 5eei = gssuj! ■ (74) 

The physical states are not u; 8 and lo\ but the ui meson and the meson which 
are mixed states 

uj = u 8 cos 9 — u 1 sin 9 

4> = Lu s sin 9 + c<Ji cos 9 . (75) 

For the vector meson it is known, that the meson is nearly a pure ss state 
as it decays mainly to kaons. Under this condition, the mixing is called ideal 
and the mixing angle is tan# = l/y/2, 9 pa 35.3°. 

If the nucleon does not couple to the (purely strange) meson, then 

9l = = s/6g NNp . (76) 

In SU(6) symmetry, which is a special case of SU(3) symmetry, only the F- 
type coupling remains for the vector mesons. This is in accordance with the 
vector dominance model [64]. In this case, the relations for the vector mesons 
can be given as in table 4. These relations reflect simple quark counting rules. 
So the A couples to the u meson only 2/3 as strong as the nucleon, as it has 
only two light quarks. On the other side, the doubly strange baryon E couples 
twice as strong to the (hidden) strange meson <fi compared to the A which has 
only one strange quark. 

The scalar coupling constants to the <j*-field are fixed by the potentials 

[4 S) « UP « 2Ui A) « 2U ( A A) « 40 MeV (77) 



27 




10 







0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

gvA/gvN 



Fig. 9. The matrix element ABaa as a function of the vector coupling constant in 
model 1 and 2. 

which is motivated from one-boson exchange models and the measured strong 
AA interaction [48] . Note that the nucleons do not couple to these new fields. 
This fixes all the parameters of model 2. 

This extended model 2 is now closer to the experimental value of the AA 
interaction matrix element (see fig. 9). Due to the additional attractive forces 
coming from the hidden strange meson exchange, the AA matrix element is 
now much stronger than in model 1. 

3. 5 Properties of Strange Hadronic Matter 

In the following we study the properties of metastable exotic multihypernu- 
clear objects (MEMO's). Metastability means that strong processes are for- 
bidden, so that the system can only decay weakly and lives therefore on the 
timescale of the weak interaction of 1CT 10 s. Inside a bound system some of 
the strong processes can be energetically forbidden due to medium effects, i.e. 
the strong reaction is Pauli-blocked. Only two processes 




A + A 
A + A 



(AE = 23MeV) 
(AE = 28MeV) 



(78) 
(79) 



have such low Q-values that this effect occurs in a bound system. 
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Fig. 10. The Q-value of the strong process AA — > EN for the system with A = 7 
and S = —4 as a function of the potential depth of the H in nuclear matter. The 
process is not allowed for negative values, hence the system can only decay weakly. 

Therefore we study systems composed of {p, n,A, S°, S~} and search for metastable 
combinations by the following recipe: 

• Start with a deeply bound nuclear core, e.g. 4 He or 56 Ni, so that the neigh- 
boring nuclei are much less bound. 

• Fill up the A-levels so that the products of the above reactions are Pauli- 
blocked. 

• Add as many EE's as possible unless the strong reactions can open. 

The effect of the Pauli-blocking is demonstrated in Fig. 10 for the lightest 
possible system H o AA 7 He. The Q-value is given by 

Q = ms + mN -2m A - B s (ls) - B N (ls) (80) 

and is plotted versus the potential depth of the S in nuclear matter. The 
system is likely to be bound and metastable if this value gets negative. This 
happens for Us > 22 MeV which is in accordance with experiment. 

Heavier systems with a doubly magic nucleon core like 208 Pb can be also filled 
up with hyperons (A, S°, and S~) selfconsistently so that they are metastable 
according to the recipe outlined. The presence of many hyperons will change 
the sequence of single-particle energies and accordingly the magic numbers 
appreciably. One gets a new shell ordering for a large strangeness \S\. Fig. 
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Fig. 11. The single particle energy of the system Pb with 70 A, 18 S , and 70 
H~ (middle column: neutron levels, right column: H° levels). The neutron levels for 
208 Pb are shown in the left column for comparison. 

11 shows the single particle energies of the system 208 Pb with 70 A, 18 S°, 
and 70 S~ which is bound by E B /A = —12.9 MeV. The neutron levels in 
the ordinary nucleus (left column) and in the strangeness-rich system (middle 
column) are compared to each other. The shell levels show some reordering 
especially for the higher lying levels due to the presence of many hyperons. For 
the protons, the change is more pronounced due to Coulomb effects and the 
magic number is now 70. The single particle levels of the S° indicates a quite 
small spin-orbit splitting. They are heavier than nucleons and in addition feel 
a smaller spin-orbit force than nucleons due to the smaller coupling constants. 
The magic numbers here are the ones of the Wood-Saxon potential without 
spin-orbit splitting: 2, 8, 18, 20, 34, 58, and 92. 

The density distribution of the system with hyperons is plotted in Fig. 12. The 
density distribution for the nucleons is still around normal nuclear density. The 
total density is increased to about 2po as the hyperon density has to be added 
up. Note, that the hyperons are treated as distinguishable point particles. 
Effects which take into account the substructure of the hadrons might alter 
this picture but will be not discussed here. The A density distributions shows a 
broad surface as the A single particle levels are occupied until the least bound 
level which ensures the Pauli blocking mechanism for the S's. The H~'s are 
less bound than the protons so that the charge density changes sign at the 
surface and gets negative. 
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Fig. 12. The density distribution of the system 2US Pb with 70 A, 18 H u , and 70 
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Fig. 13. The single particle energy of the system 180 Th with 92 A, 92 H°, and 70 
H~ (middle column: neutron levels, right column: H° levels). The neutron levels for 
180 Th are shown in the left column for comparison. 
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Fig. 14. The binding energy per baryon versus the baryon number A for Ni with 
hyperons added in model 1. 

An analogous example for model 2 with the strongly attractive hyperon- 
hyperon interaction is the system 180 Th with 92 A, 92 S°, and 70 S _ . Now the 
binding energy is Eb/A = —21 A MeV emerging from the additional attractive 
forces between the hyperons. The nucleus 180 Th itself is unstable and emits 
protons immediately. The presence of the hyperons, especially the negatively 
charged S _ 's, shift the unbound proton levels down, so that all the proton 
levels are bound. The magic numbers for nucleons has changed then from 82 
to 90 for both, neutrons and protons for large \S\. The big level changes can 
be seen in Fig. 13. In particular, the n — 1 (nodeless) levels are more bound 
than in ordinary nuclei. Also for the S levels, the nodeless li level now appears 
below the 3s level contrary to the ordinary oscillator shell ordering. 

Now we discuss the calculations for several sequences starting with a nuclear 
core and subsequently filling up the levels with hyperons. Fig. 14 shows the 
case for the nuclear core 56 Ni in model 1. When the p-shell of the A's is filled 
(Na = 8) it is energetically favourable to add S's to the system. If the next 
shell is also filled up (Na = 14), the maximum number of added S's can be 
even higher enriching the system with a lot of strangeness. 

Fig. 15 shows the sequences for the doubly magic cores 56 Ni, 132 Sn, 208 Pb, 
and 310 G (Z = 126, N = 184). The glue-like effect of the hyperons result 
in an enhanced binding energy. Adding hyperons stabilizes the systems, i.e. 
enhances the binding energy due to the opening of a new degree of freedom. 
Mass numbers of A 500 can be easily reached with a rather large binding 
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Fig. 15. The binding energy per baryon in model 1 for various sequences starting 
with 56 Ni, 132 Sn, 208 Pb and the superheavy nuclear core 310 G with 126 protons and 
184 neutrons. 

energy of Eb/A ps —13 MeV. These systems are stabilized due to the presence 
of the S~, which avoid the instability due to the Coulomb repulsion. The Pb 
sequence has systems with a total charge of Z = 12 (82 protons and 70 S~). 
For the Ni sequence, one gets also zero charged systems with 28 protons and 
28 S~. Adding further hyperons for the small systems results in less binding, 
therefore the parabolic shape for the Ni sequence. One sees a linear shape for 
the heavier cores by adding more and more hyperons. The sequences of the 
heavier systems stop when the A levels are filled up and this happens before 
there are more S~'s than protons so that one sees only a linear curve. 

The analogous sequences for model 2 are depicted in Fig. 16 for the nuclear 
cores 56 Ni and 180 Th. The other cores ( 132 Sn, 208 Pb, and 310 G) are not doubly 
magic anymore for large \S\. The system 180 Th is considerably stabilized with a 
large injection of hyperons. About 150 hyperons are needed to achieve stability. 
The region of stability extends from A = 330 to A = 578. 

Purely hyperonic systems, composites of {A, S°, H~} only, are stable in model 
2 up to A 200. They are denoted as Y which should not be confused with 
the ordinary Ytterbium nucleus. Typical binding energies are E B /A = —5 to 
—8.5 MeV. The most stable system is the one for A = 60 with 14 A, 28 S°, and 
18 S~. The Coulomb repulsion prevents systems larger than A ps 200. Inter- 
estingly, hyperonic matter is always highly negatively charged. These systems 
are lying between —0.5 < Z/A < —0.2. Moreover, the strangeness fraction is 
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Fig. 16. The binding energy per baryon in model 2 for various sequences startin 
with 56 Ni, and 180 Th. The symbol Y denotes purely hyperonic systems. 
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Fig. 17. The single particle energy of the purely hyperonic system {18A, 34H°, 28H 
in model 2. 
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Fig. 18. Binding energy versus the strangeness fraction for SHM in model 1 (nuclear 
cores 56 Ni, 132 Sn, 208 Pb, 310 G) and model 2 ( 56 Ni, 180 Th, Y denotes purely hyperonic 
systems) . 

also very high f s > 1.6 for this type of matter. For isospin saturated systems 
of {A,S°,S~} one gets the values f s = 5/3 and Z/A = —1/3, respectively. 
The lightest bound system is likely to be {2A, 2H°, 2S~} where every hyperon 
fills up their respective ls-shell. Fig. 17 shows as an example the single par- 
ticle energy for the hyperonic system {18A, 342°, 28H~}. The ls-state has a 
binding energy between 50 and 60 MeV for all three hyperon species. This 
reflects more or less the assumed universal hyperon-hyperon interaction. The 
energy levels of the S° are a little bit deeper than the ones of the A due to 
the higher mass. The Coulomb repulsion shifts the states of the S~ to smaller 
binding energy compared to the others. 

The binding energy of strange hadronic matter is plotted in Fig. 18 versus 
the strangeness fraction f s = \S\/A for both models for various nuclear cores. 
It is interesting to note that the minima of the curves are located around 
f s ~ 0.6 for model 1 and f s ~ 1.0 for model 2 which is in the same range 
as for strangelets. The minima found in Fig. 18 correspond to charge to mass 
ratios around zero very similar to the property of strange quark matter in 
its groundstate. Nevertheless, the main difference between strangelets and 
strange hadronic matter, except, of course, of their internal structure, is that 
the maximum binding energy of Eb/A r* —13 MeV (model 1) and Eb/A r* 
—23 MeV (model 2) is by far to less to overcome the mass difference between 
nucleons and hyperons of about 177 MeV. Strange hadronic matter will thus 
decay weakly on a timescale comparable to the lifetime of a A-hyperon. 
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4 STRANGE MATTER IN RELATIVISTIC HEAVY ION COL- 
LISIONS 



Collision experiments of heavy ions at high bombarding energies may be the 
only unrevealing chance and opportunity for both the production and detec- 
tion of small pieces of strange quark matter: One major goal of these exper- 
iments is to unfold the temporary creation and existence of the quark gluon 
plasma phase. A direct proof of this fleetingly small moment turns out to 
be a difficult and complicated task. One therefore has to find signals which 
come directly and only from the transient existence of the QGP. This will 
be different, if a (met a-) stable strange quark droplet will show up during the 
break-up. This is quite similar to the situation in the early universe [11], where 
a remnant (or 'ash') of the QGP is created. In fact, it seems rather selfevident 
that a deconfined quark matter state could probably act as a good starting 
point for the agglomeration of quarks to strangelets, because the quarks are 
thought to move quasifreely over longer distances similar as in the situation 
inside strange quark matter. 

In this section, we review the ideas that had led to the search for strangelets 
at present heavy ion experiments at Brookhaven and at CERN: In the first 
subsection, the important mechanism of strangeness separation in the phase 
transition of a quark gluon plasma back to hadronic degrees of freedom will 
be explained [3]. The second subsection then addresses the exciting possi- 
bility of producing strange quark matter droplets in heavy ion collisions. In 
particular we will discuss why such droplets might cool and form long-lived 
cold strangelets. In the last subsection we will close our discussion in criti- 
cally emphasizing the detection possibilities of strangelets and small MEMOs 
by their properties and lifetimes, also in respect to the present experimental 
undertaking at Brookhaven and at CERN. 



4-1 Strangeness separation during the QGP phase transition 

In nuclear collisions, strangeness can only be produced in ss pairs due to the 
conservation of hypercharge in strong interactions. Thus, at first sight, there 
seems to be no chance to really succeed for producing strangelets: The net 
strangeness of the plasma state counting the difference of strange to antis- 
trange quarks is zero from the onset of the fireball's expansion. In this sub- 
section we want to give arguments, that first if a baryon-rich and hot QGP is 
created in such collisions and second if the strangeness degree of freedom is 
nearly saturated, the strange quarks will separate from the antistrange quarks 
in a nearly equilibrium deconfinement <-> confinement transition and will pre- 
dominantly remain in the plasma phase [3]. 
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According to fluid dynamic or microscopic models baryon stopping, (baryon- 
) dense and hot matter formation and partial thermalization has to be expected 
at AGS and CERN energies for central collisions and heavy systems [65] . One 
should be aware of the fact, that the energy and baryon densities reached are 
so large, that a phase transition to a quark gluon plasma may have occurred. 

It was suggested and shown already some time ago that an abundant number of 
strange and antistrange quarks are produced in a hot QGP [66] by gluon fusion 
and that, accordingly, strangeness saturates the phase space after a very short 
equilibration time which may actually be shorter than the lifetime of the QGP 
phase [6,67]. In the meantime, an enhancement of strange particle production 
in nuclear collisions has been observed in many experiments [68-71]. However, 
it has also been learned that such an enhancement alone does not make a 
reliable signature for the QGP. Strange particles, especially K mesons and A 
hyperons can be copiously produced in hadronic reactions before the nuclear 
fireball reaches equilibrium [72]. Yet, the enhancement of the A hyperon over 
a wider rapidity range [71] and especially the observed and strongly enhanced 
yields of antihyperons [69-71], which were also proposed as a signature for 
QGP [6], require particular microscopic modifications of the hadronic cascades: 
Color rope formation [73]; multiple string breaking and (!) decaying multi- 
quark droplets [74]. Such modifications can already be considered as precursor 
phenomena associated with QGP formation. 

In the following, the deconfinement phase transition is assumed to be of first 
order, implying that the relaxation times for chemical transmutations as well 
as the hadronization time are small compared to the overall transition time. 
The Gibbs criteria read 



Tqgp 


— Thg 


Pqgp 


= Phg 


b 

Vqgp 


B 

— f^HG 


Vqgp 


= ^HG 



In particular the last condition is significant for the separation to be discussed: 
In a Gibbs phase equilibrium the chemical potentials are continuous across the 
phase boundary, whereas the corresponding densities are (or might be) dis- 
continuous. This first implies that the abundance of especially the strange 
hadrons during the equilibrium transition are governed by the same quark- 
chemical potentials, ji q for the light quarks and fj, s for the strange quarks, 
e.g. 

^K,K — "F Us , 
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fi E ^ = ±fi q ±2fx s , (82) 

Second, the corresponding density to the strange quark potential is the net 
strangeness content, i.e. the strangeness ratio f s defined in (3). In general 
this ratios need not to be the same during transition: f® GP (fi q ,fi s ,T) ^ 



Fig. 19. Schematic picture of the separation of strangeness during phase transition 
in a baryon-rich system 

Consider now the phase transition of the QGP to the hadron gas at some 
critical temperature. How does the strange and antistrange quarks hadronize 
during this transition? There is no stringent reason why these different quarks 
do hadronize in the same manner and time, especially if one thinks of a baryon 
rich system. At the beginning of the hadronization (t — ts) there is only the 
QGP phase. Because of an overall strangeness conservation the number of 
strange quarks equals the number of antistrange quarks, the strange chemical 
potential fi s = and the net strangeness content of the plasma phase is that 
of the system, ff GP {t = ts) = 0. Now, quarks combine to hadronic particles 
and leave the plasma phase, but hadrons may also decompose and go back to 
the plasma state. Because of the progressing expansion of the whole fireball, 
however, the hadron phase turns bigger while the plasma phase decreases. 
It is 'simple' for the antistrange quarks to materialize into kaons, because of 
the (immense) surplus of massless quarks compared to their antiquarks. The 
strange quarks can combine to A-particles, but these are rather heavy and it 
is energetically much easier for them to stay in the plasma, when hadroniza- 
tion proceeds. The possible strangeness content in the two different phases is 
schematically drawn in Fig. 19. A large antistrangeness builds up in the hadron 
matter while the QGP retains a large net strangeness excess. This separation 
will occur only when the system carries a positive net baryon number. 
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Fig. 20. Fraction of net strangeness to baryon number present in the QGP phase as 
a function of the baryochemical potential fj, q and the volume fraction \. Note that 
f s can exceed 0.5. The path of an isentropic expansion is also shown. 

To be more specific, requiring that the total strangeness during the phase tran- 
sition in the combined system of two phases vanishes, leads to the constraint 

V QG p(ps - ps) + V H g(pk + Py - Pk ~ Py) = . (83) 



Together with the Gibbs criteria of pressure equilibrium, (81), the phase tran- 
sition temperature and the strange quark potential are then implicit functions 
of the quark chemical potential p q , and, for example, the volume fraction \ °f 
the hadron phase to the total system, i.e. x = V HG /(V HG + V QGP ) [3]. In Fig. 
20 the fraction f® GP (p q , p s ,T) of strange quarks to all quarks present in the 
QGP is shown as a function of the quark potential p q and the volume fraction 
X- At the beginning of the phase transition (V HG = 0) the well known result 
p s = is recovered, s and s quarks are produced in pairs only. On the other 
hand, for vanishing QGP (V QGP = 0) zero net strangeness leads to a nonzero 
value of the strange chemical potential p s > [3]. A different strange particle 
production shows up in the dominance of the associated production over the 
direct pair production at finite baryon density. During the coexistence of the 
two phases, an additional channel opens up for the strangeness: besides the as- 
sociated pair production in the hadron gas it is possible to have, for example, 
associated production of a K - meson in the sector of the hadron gas phase 
and the s-quark staying in the QGP. This leads to a net strangeness content 
ff GP (t > tg) larger than in the quark phase and to a diminished hyperon 
abundance, f^ G < 0, the hyperons being too massive, in the hadronic gas sec- 
tor at finite net baryon densities. At vanishing overall net baryon density, i.e. 
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at n q = 0, (and vanishing net strangeness density) the transition is completely 
symmetric for strange and antistrange particles, i.e. fi s = throughout the 
transition, and separation can not occur as one would expect. 

Of special importance is the fact that the accumulation of s quarks in the 
plasma phase grows with decreasing plasma volume. This opens up the possi- 
bility that s-quarks may be bound not only in hyperons and strange mesons: 
They could form strange quark matter clusters or MEMOs which might be 
metastable objects. A baryon-rich environment may thus be more feasible for 
strangelets than the situation in the hot quark phase in the early universe with 
only a tiny small surplus. Also the cooling of the plasma should in addition 
not only be due to the expansion, but also due to 'prefreezeout' evaporation 
processes of hadronic particles, because of the relatively small length scale in 
heavy ion collisions compared to the astrophysical picture: Surface processes 
like evaporation should become important. This issue and the possible distil- 
lation of strangelets will be discussed in the next subsection. 

Finally we remark that the separation mechanism can be probed by density 
interferometry with hyperons [75] or kaons [76]. The hadrons with negative 
strangeness, the A, K and K°, are expected to be produced mainly at the 
last stage of the phase transition when the size of the quark phase volume 
has become quite small. As the hyperons or antikaons are quite heavy, their 
thermal velocity is relatively small so that at freeze-out they probably still will 
be localized near the center of the system. Such a close localization r ps 2 — 3 
fm can be seen by measuring two particle correlations with the Hanbury- 
Brown-Twiss method. In Fig. 21 the pair correlation function of two A particles 
is shown under the assumption that a rather strong possible resonance channel 
contributes at low momenta. This might occur because of a possible strong 
attractive A-A-interaction (see section 3). Especially with such a channel small 
source sizes are relatively clean to determine. (Of course, the observation of 
such a resonance would give some exciting new experimental insight on the 
interaction among two As.) If the separation does not occur or if there is no 
QGP the deduced source size should be nearly the same as for kaons or pions 
and thus considerably larger. 



4-2 Strangelet distillation 

Besides the onset of strangeness separation during the equilibrium phase tran- 
sition there exist another, somehow related argument which could lead to 
strangelet formation in ultrarelativistic heavy ion collisions [4,77]. Rapid kaon 
emission from the 'surface' (or, in other words, the outer region) of the fireball 
can result in an even stronger enhancement of the s-quark abundance in the 
quark phase. As important, this prompt kaon (and, of course, pion) emission 
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Fig. 21. The pair correlation function of emitted A hyperons is shown for different 
source sizes under the assumption of a relatively strong resonance channel. 

may in addition cool the quark phase, so that there really might be a chance 
that the quark plasma droplet condenses into metastable or stable, rather cold 
droplets of strange quark matter [4]. 

In order to model the evolution of an initially hot fireball a two phase equi- 
librium description between the hadron gas and the QGP was combined with 
the nonequilibrium radiation by incorporating the rapid freeze-out of hadrons 
from the hadron phase surrounding the QGP droplet during the phase tran- 
sition [5]. In particular we will then address the question whether, how and 
why a hot QGP droplet can cool and will then form (or better 'become') a 
long-lived cold lump of strange quark matter. Before we turn to the model we 
summarize the reasoning for the existence and the possible consequences of 
the evaporation processes. 

Let us make some crude estimates about the qualitative features of the surface 
radiation of a baryon-rich fireball. One first might think of meson radiation off 
the initial pure quark phase by various microscopic processes [77-80]. The most 
dominant particles to be evaporated are the pions. For the strange mesons the 
K + , K° are more easily radiated in a baryon-rich environment off the hot 
surface of the quark phase than K~ and K° [4]. Similarly and probably even 
more important should be the meson radiation while the system is in the phase 
coexistence region. On one hand most of the total collision time is spent in 
this region On the other hand, the possibility of reabsorption of mesons is 
very unlikely, since the system approaches freeze-out. Meson evaporation in 
both stages of the expansion just described carries away entropy, energy and 
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antistrangeness. Therefore, the residual expanding fireball, which is in the 
mixed phase, loses entropy and is charged up with a net strangeness larger than 
zero. For a rough estimate a thermal black-body-type radiation formula can 
be employed. Accordingly the emission of Ks containing an g-quark are then 
suppressed from baryon-rich matter by a factor ~ e ~ 2fMq / T as compared to the 
Ks. Adding the resulting numbers over the different stages of the fireball, a loss 
in entropy per baryon of about 15 units [4] and a net strangeness enrichment of 
the total system of about fl° tal = 0.75 is found [77,4]! The assumption of black 
body radiation seems intuitively plausible from detailed balance arguments 
[81] if one thinks that each hadron moving towards the 'surface' of the plasma 
phase will first coalescence and then be completely absorbed. 

The mechanism by which a QGP state is converted into hadrons is a major 
uncertainty in the different descriptions. The hadronization transition has of- 
ten been described by geometric and statistical models, where the matter is 
assumed to be in partial or complete equilibrium during the whole expansion 
phase. A hydrodynamic expansion is often assumed, which is substantially 
modified by the phase transition [82,67]. Partial departure from (strange) 
chemical equilibrium is allowed in the flavour kinetic model [6,67,83]. These 
rate calculations suggested that chemical equilibrium is reached in the hadron 
phase, if the system evolves from the deconfined phase: gluon fusion yields 
fast equilibration of strangeness. The subsequent re-hadronization affects that 
the strange hadron yields closely approach the hadron equilibrium yields from 
above. 

A more realistic scenario must take into account the competition of the col- 
lective expansion with the particle radiation from the surface of the hadronic 
fireball before 'freeze out'. In a second class of models [84,85] it is assumed 
that the plasma breaks up into droplets at some point during the hadroniza- 
tion transition. In a 'cascade' model [85] the emission and reabsorption of 
pions from the surface of the hot plasma was studied. For the present pur- 
pose, evaporation of other hadrons must also be taken into account. Nucleons 
and strange particles (e.g. kaons and hyperons) are particularly important. 
One important outcome of the cascade model [85] was the observation that 
the total entropy is approximately constant during the hadronization. 

Our following model [5] combines in a simple, yet plausible way the two dif- 
ferent classes. The expansion of the QGP droplet during the phase transition 
will be described within a two-phase equilibrium; in particular the strangeness 
degree of freedom is thought to stay in chemical equilibrium because the com- 
plete hadronic particle production is driven by the plasma phase, as suggested 
by the rate calculations. The competition of the collective expansion with 
the particle radiation is incorporated by rapid freeze-out of hadrons from the 
outer layer of the hadron phase surrounding the QGP droplet at the Gibbs 
phase transition point. The expansion is here assumed to be quasiisentrope, 



42 



K 




7T 



Fig. 22. Hadron fluid surrounds the QGP at the phase transition. Particles evapo- 
rate from the hadronic region. New hadrons emerge out of the plasma by hadroniza- 
tion. 

as suggested by the cascade model. This scenario is visualized in Fig. 22. 

The global properties of the two phase system then change in time in accord 
with the following differential equations for the baryon number, the entropy 
and the net strangeness number of the total system: 



—A tot = — T A HG 
dt 

^S tot = -TS HG (84) 
at 

where T = ( j is the effective ('unique') rate of particle (of con- 
verted hadron gas volume) evaporated from the hadron phase. These three 
equations may easily be combined in a physically illustrative form: 



A / C\tot 

J =-TX HG (l-X HG 



dt \A 

d 
~dt 




/r = -r \ HG (i - \ HG ) (/f G - ff GP ) , (85) 



where \ HG = 1 — \Q GP = A HG /A tot is the ratio of baryon number contained in 
the hadron phase over the total baryon number in the system, f s defines the 
total net strangeness content of the system, and 4 is the entropy per baryon. 
The three equations (84,85) constitute a set of coupled differential equations 



43 



for A tot , (S/A) tot and /*°*. The latter two are given in terms of the hadron and 
quark phase content by 

, Ox tot , Q X QGP , qxHG 

( I ) = (^J (/V ^ T) (1 - A" G ) + (-J ,u s , T)\ HG (86) 

/r = /f GP (^, ^, T) (1 - A HG ) + f™(» q , ^ T)X HG . 

The equation of state used consists of the Bag-model for the QGP and a mix- 
ture of relativistic Bose-Einstein and Fermi-Dirac gases of the well established 
nonstrange and strange hadrons up to M ~ 2 GeV for the hadron matter. 

Although the baryon number A and the strangeness f s are conserved under 
strong interactions, their value in the system changes with time according to 
the above equations due to the evaporation process. This calculation requires 
solving simultaneously the equations of motion and the Gibbs phase equilib- 
rium conditions (81), which specify the intrinsic variables (e.g. the chemical 

potentials) of the functions (f s j QGP/Had ^ T ) anc i fQ GP / Had ((j, q , fi s ,T). 
Still the effective rate T and the baryon fraction \ HG are not specified. Intu- 
itively, the plasma is surrounded by the hadron phase, thus we take it as a 
shell of constant thickness of Ar ~ 1 fm (this is an ad hoc assumption, how- 
ever, as it turns out, the results depend only weakly on this thickness). This 
then specifies X HG . Some unique hadron gas volume of the outer layer is then 
evaporated in a small time interval At. For the unknown evaporation time we 
use an averaged rate of all particles at the considered temperature, dictating 
the time-scale for the evolution of the hadronic particles out of the hadron 
phase. Altogether, the quasiadiabatic expansion during hadronization and the 
evaporation of the hadrons in the most outer regions are both incorporated 
and dictate the evolution of the (plasma) droplet and its intrinsic variables 
fi q ,fi s and T. 

For the following arguments it is important to study the expansion of an ini- 
tially hot QGP fireball with special attention to the evolution of the strangeness 
during the phase transition. As outlined above, early non-equilibrium particle 
radiation off the pure hot and dense (baryon-rich) QGP fireball is expected to 
be important in the initial phase of the reaction. Pions and kaons (containing 
an s-quark) (and only a minor number of antikaons) are emitted from the 
surface of the plasma. Thus the QGP is enriched with finite net strangeness 
even before the phase transition point is reached. The net strangeness enrich- 
ment resulting from the early black body radiation off the pure QGP-droplet 
has been estimated to be in the range 

jinit < 5 [ 77j4 j_ Thege valueg 

nave 

been used in the present model as input for the initial condition f s (to) when 
hadronization starts. In a complete isentropic expansion with a total (initial) 
net strangeness content but without further particle evaporation it was (al- 
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ready) shown [4,86] that for not so high initial entropies the evolution can end 
with a small and (meta-)stable strangelet, if the employed QGP equation of 
state will allow for such exotic configurations to exist at small temperatures. 
We will see in the following how things will work. 



Fig. 23. (a) Baryon number, strangeness content and temperature of the quark glob 
during complete hadronization as a function of time for a very large bag constant 
B~* = 235 MeV. The initial values are an initial baryon content of Ab(£o) = 100, 
an entropy per baryon ratio of ^(to) = 25 and an initial net strangeness fraction of 
fs(to) = 0.25. Note the strong increase of the strangeness content with time, (b) The 
same situation as in (a), however, for a small bag constant B± = 145 MeV, when a 
strangelet is distilled. One observes a strong decrease in the evolving temperature. 

Fig. 23 gives an impression how the hadronization proceeds for a large bag 
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constant (.B 1 / 4 = 235 MeV - no strangelet in the groundstate) and a small bag 
constant (i? 1 / 4 = 145 MeV). The initial parameters are a small net strangeness 
content of f s (to) = 0.25 and a moderate entropy per baryon ratio of -f(io) = 25 
(which is expected at CERN SPS energies): For the large bag constant the 
system hadronizes completely in t ~ 8 ^p, which is customarily expected 
and thus not surprising. The quark droplet remains unstable until the strange 
quarks have clustered into A-particles and other strange hadrons to carry 
away the strangeness and the plasma has completely vanished into standard 
particles. Yet, a strong increase of the net strangeness of the system is found in 
both situations, which is basically a confirmation of the strangeness separation 
mechanism and will be analyzed in some more detail below. The plasma drop 
reaches a strangeness fraction of f s > 1.5 when the volume becomes small. 
Indeed, for the small bag constant, however, a cold strangelet emerges from 
the expansion and evaporation process with an approximate baryon number 
of A ~ 22, a radius of R ~ 2.5 fm, and a net strangeness fraction of f s ~ 1.5, 
i.e. a charge to baryon ratio Z/A ~ —0.25! This would be an interesting object 
also from its global properties: It would comprise a nucleus of positive baryon 
number, but negative charge. 
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Fig. 24. The number of emitted particles are shown versus time in a situation, 
when a strangelet is distilled. The conditions are the same as in Fig. 23, however, 
no initial net strangeness is assumed. 

The particle yields of the most frequently emitted species are shown in Fig. 24 
for the case, where a strangelet condenses after ~ 100 — , although the total 
net strangeness was zero at the beginning of the hadronization. Their time 
dependence is more or less as expected: Roughly 300 pions and 32 kaons are 
produced prior to most of the other hadron species, followed by the nucleons 
(~ 77) and at the late stage by some antikaons (~ 3) and hyperons (~ 6). 
The total K + / 7r + -ratio is then approximately 0.16. The surplus of 17 in baryon 
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number and of 23 in strangeness are contained in the missing strangelet, which 
has to have a strangeness fraction of ~ 1.4! The temporal behavior of the 
emitted particles underlines the important role played by the kaons (and pions) 
during hadronization. The early leaving kaons do enrich the remaining system 
with net strangeness. The pions carry about half of the available entropy off 
the system. For this reason one might expect that the nonequilibrium aspects 
can become even more efficient at still higher entropies, when more kaons and 
pions are produced. 

One observes that in both cases depicted in Fig. 23 the net strangeness in- 
creases. However, in the first case, at high critical temperatures, the tempera- 
ture of the system increases slightly, although one might have expected some 
cooling by the emission of particles, whereas in second case, strangelet forma- 
tion goes hand in hand with strong cooling. As a compelling consequence the 
temperature (and entropy) decreases drastically throughout the hadronization 
and condensation process. Thus the pressure goes to zero and an absolutely 
stable, cold strangelet is created. The energy per baryon of the remaining 
strange quark droplet converges to its zero temperature value. 

Why does in both (model) cases the net strangeness content increase? Consider 
therefore the physical content of eqs. (85). Obviously, T, \ HG and (1 — X HG ) 
are larger than zero. Therefore, the total net strangeness fraction /*°* of the 
system will always increase or decrease, respectively, with time, if the frac- 
tion f^ G is smaller or larger, respectively, than f® GP . Of course, this is just 
the question whether or not the separation of strangeness does occur. The 
abundant presence of the kaons enforces f^ G < early on! Also later on f^ G 
remains smaller than f® GP , and therefore /*°* increases as function of time 

d/r>o). 

Why is there cooling in the one and reheating in the other case? Here a rather 
similar argument holds for the total entropy. The total entropy of the remain- 
ing system (S/A) tot will decrease (4(5 / A) tot < 0), i.e. the system will be 
cooled, if and only if the specific entropy per baryon in the hadron phase ex- 
ceeds that in the quark phase, (S/A) HG > (S/A) QGP . The system will reheat, 
i.e. (S/A) tot will increase (f t (S/A) tot > 0), if and only if {S/A) HG < (S/A)^ GP . 
It was pointed out already in ref. [4] that the strangelet formation can only go 
hand in hand with strong cooling rather than reheating. This is the case when 
the bag constant B± is small and allows for the existence of (meta-)stable 
strange quark matter states. Although this intimate connection between cool- 
ing of the plasma phase and the existence of strange quark matter is intriguing, 
it might be valid only within the simple parametrization of the quark gluon 
plasma phase within a bag model description. Ultimately, whether (S/A) HG 
is larger or smaller than (S/Ap at finite, nonvanishing chemical potentials 
might theoretically only be proven rigorously by lattice gauge calculations 
in the future, as also the principle existence of (meta-)stable strange quark 
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matter. 



In the situation when a strangelet gets distilled out of the original QGP phase, 
one observes [5] that the values of the two chemical potentials approach 'rather' 
quickly their asymptotic numbers, even when the temperature of the strangelet 
is still not so small, T ~ 40 — 60 MeV: or metastable): 

• fi q ~ 200 - 275 MeV , fi s ~ 300 — 325 MeV (B 3 = 145 MeV); 

• fi q ~ 200 - 250 MeV , fi s ~ 350 — 400 MeV (B* = 160 MeV). 

The evaporation of the particles forces these potentials to become quasista- 
tionary; a detailed balance for the emission of nucleons and antikaons has 
adjusted itself. The hyperon emission is still strongly suppressed compared to 
the nucleons, although the strange chemical potential is already larger than 
the light quark potential. (For higher bag parameters, the hyperons become 
more and more important. It is exactly when strangelets are getting barely 
stable to unstable when the hyperon emission becomes as or even more effi- 
cient in the late stages of the plasma phase than the one for the nucleons.) 
A 'simple' but basic idea for the possible formation of strange quark matter 
droplets is due to behavior of the nucleons. As discussed in the section 2, the 
total energy per baryon (at finite temperature this has to be substituted by 
the free energy) can easily be lowered by assembling the non-strange quarks 
into pure nucleonic degrees of freedom as long as the strangelet is not at its 
minimum value in f s . This happens rather fast. For finite temperatures the 
surplus in thermal energy provides additional nucleon evaporation, so that 
this minimum point will be 'overshot' and the strangelet turns out to have 
a strangeness fraction larger than 1. To further illustrate how then stable 
strangelets cool and survive, we note first that with the above numbers an 
effective nucleon binding energy I n = m n — 3fi q ~ 120 — 350 MeV is found in 
the late stages of the process. The energy per baryon of a strangelet at small 
temperature and zero pressure is approximately given by 



Any late emission of a nucleon will now decrease the energy per baryon of the 



In addition, the differential of eq. (87) is A(E/A) = 27TAT. Combining these 
two simplified expressions and integrating from the initial temperature T { to 
T = yields roughly the final baryon number of the strangelet: 




(87) 



blob: 





(88) 
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For I n ~ 250 MeV and % ~ 20 MeV, a total of ~ 92 percent of the initial 
baryon number would remain in the droplet. Hence, evaporation of baryons 
from quark matter will be suppressed at very low temperatures (and will cool 
the droplet), because there is simply not enough thermal energy available to 
power the complete evaporation. The strangelet becomes cold and stable. 

We have so far demonstrated that under certain circumstances, the hadroniza- 
tion can result in the distillation of a rather cold strange quark matter droplet. 
This debris of the quark gluon plasma can, of course, only survive, if the energy 
per baryon of the strange quark object is smaller than the energy per baryon 
of the corresponding hadronic system at the same f s . For the bag model the 
strangelet is stable for = 145 MeV or metastable for B* = 160 — 180 
MeV. For realistic parameters modeling the initial situation of the QGP fire- 
ball either in the stable or metastable situation one might expect a reasonable 
strangelet with a baryon number ~ 10 - 30 [5]. Negatively charged strangelets 
result for all cases calculated in [5]: 

• /,(*-> W) ~ 1.1 — 1.5, Z/A ~ (-0.05) - (-0.25) - stable case ; 

• f s (t -> 'oo') ~ 1.5 - 2.0 , Z/A ~ (-0.25) - (-0.5) - metastable case . 

One should emphasize the fact that the charge of the distilled strangelet (with 
positive baryon number) is found to be negative! This seems counterintuitive 
as the most stable configurations should actually be, though small, positively 
charged. This has to do with the final nucleon emission which is still energet- 
ically possible if the strangelet carries still a positive charge. Its charge then 
can only be changed by weak processes. (Please note our previous discussion 
above. We will come back to this point also in the next section when we ad- 
dress the important question of what are the most likeliest candidates to be 
seen in the dedicated experiments.) 

As long as finite size effects like surface tension or curvature contributions 
are neglected, the equations of motion (84) are scale invariant. All extensive 
properties, the rate T and the time t do scale for example with the initial 
baryon content A B (t ). Thus, in order to distill a strangelet with a baryon 
number much larger than 1, the initial quark gluon plasma droplet must be 
fairly large. (Not only) for this reason the heavy ion experiments using the 
heaviest ions, like the available Au- and Pb-beams, should be most favourable 
for strangelet search. 

The distillery works even for larger initial entropies S/A=50 or 100 [5,87]. A 
high initial entropy does not necessarily prohibit strangelet formation. Abun- 
dant kaon production enriches the plasma rapidly with net strangeness at 
high entropies. This might offer to look for strangelet production at the high- 
est bombarding energies available in the present and future for very heavy 
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systems, e.g. at the CERN SPS (E LAB ~ 200^-) or even at RHIC (E LAB 
~ facilities. One might even think that the distillation indeed could 

also work for very small initial chemical potentials, i.e. \i q jT <C 1, as will be 
expected in the central rapidity region at RHIC energies (although there still 
might be a noticeable, nonvanishing net baryon excess) and LHC energies [87]. 
There might be fluctuations at some, but no particular rapidity region where 
a small but net excess in baryon number as well as in strangeness number 
statistically occurs. Within the present model this can lead to a distillation of 
very small strangelets of a size A B < 5 [87]. There is a lot of wishful thinking 
in this scenario, as it assumes that at the onset of the transition the system 
is in chemical equilibrium and due to the fluctuations is described by very 
small, but finite chemical potentials. It is amusing to note some analogy with 
the cosmological scenario proposed by Witten [11]. The tiny net baryon excess 
in the early universe hiding in the 'high-temperature' phase is tremendously 
small, even compared to the numbers expected at the LHC energies due to 
fluctuations. The distillation of strange quark matter Witten visualizes as a 
shrinking, leaking 'balloon', where only neutrinos are allowed to escape, leav- 
ing the net baryon density back and trapped inside the balloon. This picture 
assumes neutrino losses are the main way for the high-temperature phase to 
lose energy and baryon diffusion to be negligible, while in fact neutrino losses 
and surface evaporation might appear comparable. In the present case the 
mesons take the part of the neutrinos and it is the question how efficient their 
emission is compared to the evaporation of baryons, so that a cooling really 
might take place. 

We conclude that there are two essential processes which can favour the forma- 
tion of strangelets from a baryon-rich QGP formed in ultrarelativistic heavy 
ion collisions. The first (as a prerequisite) is the s-s-separation mechanism 
discussed in the section before. This leaves us with a quark phase in the 
coexistence region of hadronic and quark matter, which is charged up with 
strangeness. The s-quarks, created in pairs with s-quarks in the early quark- 
gluon fireball, remain in the quark phase during the phase transition, the 
s-quarks materialize mainly into kaons. The second is the evaporation of the 
hadronic gas with its antistrangeness excess, which charges the remaining sys- 
tem with net strangeness. Besides the expansion additional pion and nucleon 
evaporation should help to allow for a possible, yet necessary cooling. All these 
processes are statistical in nature. These calculations are strictly valid only 
for large systems. Large fluctuations around these average predictions have to 
expected, in particular for the strangeness production per event and for the 
evaporation processes. In a real heavy ion collision an idealized situation as put 
forward in the present description is probably never reached. The results here 
should be seen with this in mind. They should illustrate that if (meta-)stable 
strangelets do (or can) exist in nature, there might really be a chance for their 
production at the present and future relativistic heavy ion experiments. 
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4-3 Detectability in heavy ion experiments 

An important prediction of the exploited model in the previous subsection is 
the negative charge to mass ratio, Z/A ~ —0.1 for absolutely stable strangelets, 
and Z/A ~ —0.45 for metastable droplets. Still, not restricting too much 
to this model, the produced strangelets or MEMOs might also be slightly 
positive [4]. It then could be detected by its unusual charge-to-mass ratio 
(~ +0.15 < Z/A < -0.5). 

On the other hand, simple coalescence estimates give production probabilities 
of strange clusters of the order of 10 3_As_ l s l, where S denotes the strangeness 
and A B the baryon number of the cluster [9,88]. Small clusters with A B + \S\ < 
r + 3, where r is the sensitivity of the apparatus (presently r < 12), are most 
favoured for detection. Therefore, if strangelets or MEMOs are formed due to 
this scenario, baryon numbers of Ab < 12 are expected. 

It is important to note that these objects are a new form of matter, not a 
specific new particle. The strange droplets produced in these reactions do 
not come in the form of a single type of particle. Many different sizes of 
droplets may be produced, spanning a range in mass, charge, and strangeness 
content. The experimental task of finding the new form of matter is therefore 
challenging. Here any detected particle having an unusual charge to mass ratio 
is a potential hypermatter candidate. 

To identify a particle or cluster, its charge and mass need to be measured. To 
determine that the particle is a new form of strange matter, its strangeness 
content must also be revealed. The experimental approach is first to find 'ob- 
jects' having a peculiar or new charge/mass ratio. (The strangeness might 
be seen by interaction with a secondary nucleus : multiple production of As, 
Ss, Ss and Ks in such a secondary reaction would signal its existence.) The 
key idea here is that the charge/mass ratio will be unlike that of any normal 
nuclear isotope (the 8 He with a Z/A = 0.25 would be the isotope candi- 
date with the smallest ratio). Strangelets (or MEMOs) would have a charge 
~ 0, being slightly positively or negatively charged. In particular in the range 
—0.25 < Z/A < +0.25 there exists no quasistable form of nuclei or antinuclei. 
Such a range will has been covered by the E864[89] experiment at Brookhaven. 
E878, the successor of E858[90], using a focusing spectrometer at zero degree, 
is seizing a much smaller selected range, which in respect to cover still the full 
range of interest can be steadily adjusted. A similar technique like in E878 is 
also been employed by the Newmass collaboration[91] (NA52 experiment) at 
CERN at much higher energies, and at the present time also with the heavy 
P6-beam. 

The question remains, on which time-scale weak decay or flavour changing 
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modes will appear. How would a strangelet or MEMO then look like when 
passing through the detector? Employing TOF-techniques to reveal the veloc- 
ity and thus the charge to mass ratio, the experimental setup sets a natural 
time scale ~ 1CT 8 sec. So, an important question we finally have to address 
are the lifetimes of these objects. In the following, we subsequently discuss 
the properties of both forms of strange matter and the possible long- and 
short-lived candidates referring to [92]. 

Whether or not strangelets exist depends crucially on the value of the bag 
constant which is not known for such strange and big systems. For a bag 
constant of B 1 ^ = 145 MeV, the original value of the MIT bag model fit, 
strangelets are absolutely stable, for bag constants up to B 1 ^ m 180 MeV 
strangelets are metastable, i.e. they can decay by weak interactions. So any- 
thing between absolutely stable and unbound is possible. Nevertheless, for the 
following arguments one needs only three basic assumptions: 

(1) Strange quark matter is at least metastable. 

(2) There exists a local minimum for the total energy per baryon of strange 
quark matter at a finite strangeness fraction f s = \S\/A. 

(3) The relativistic shell model can be used for strangelets. 

With these assumptions we demonstrate that there exists a valley of stability 
at low mass numbers and that these strangelets are highly negatively charged 
contrary to former findings. 

The MIT bag model is used here as a guideline only. Fig. 25 shows the energy 
per baryon number of isospin symmetric strangelets as a function of f s for 
A < 40 for a bag parameter of _B 1//4 = 170 MeV. Now there are three different 
processes which will shift a strangelet emerging from a heavy ion collisions to 
a very high strangeness fraction. First, the strangelets sitting above the line 
drawn between the nucleon and the hyperon masses will decay to a mixture 
of nucleons and hyperons by strong interactions completely as this is ener- 
getically favored. Second, the strangelets located between that line and the 
tangent construction starting at the nucleon mass (denoted as strong) can 
still decay strongly on a timescale of 10~ 20 sec by emitting nucleons and hy- 
perons. They will be shifted to a higher strangeness fraction until they reach 
the tangent point at f s ~ 1.4 (confer our discussion at the end of the previous 
subsection). Third, weak nucleon decay can occur for the strangelets between 
the former tangent and the other tangent (denoted as weak) starting at the 
nucleon mass and f s — 1 (as weak interaction change one unit of strangeness) 
[9]. The timescale for a weak nucleon decay has been estimated to lie between 
10~ 7 -10~ 9 sec [93]. More conservatively, one could also argue that it might 
be comparable to that of a A-particle. For a strangelet with f 8 >l the weak 
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Fig. 25. The energy per baryon E/Ab of isospin symmetric strangelets with Ab < 40 
for a bag constant of B l / A = 170 MeV versus the strangeness fraction f s . The solid 
line connects the masses of nucleon, A, S and and stands for free baryon matter. 

nucleon decay will enhance the strangeness fraction as 



Hence, strangelets surviving strong and weak nucleon decay (and which thus 
can be detected with present experimental setups) can be sitting at a very 
high strangeness fraction of f s 2.2 which is the weak tangent point in Fig. 
25. For isospin symmetric systems, this large strangeness fraction corresponds 
to a charge fraction of 

j = \ (1 - /.) = "0.55 (90) 



which indicates highly charged strange quark matter. This is contrary to the 
conventional picture that strangelets have a slightly positive charge-to-mass 
ratio which is the case for strange matter sitting in the minimum of the curve 
plotted in Fig. 25. But as pointed out before, the combined effect of strong and 
weak hadronic decay will shift strangelets emerging from a heavy ion collision 
to much higher values of f s and therefore to highly negatively charged objects! 
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This simplified picture is only valid in bulk matter. For finite systems, which 
we are interested in, shell effects will be important (see already our discussion 
in section 2). Already in Fig. 25 one sees that shell effects are at the order of 
100 MeV per baryon number! Hence, we expect that strangelets with a closed 
shell can be very deeply bound. These 'magic' strangelets are most likely to be 
stable against strong and weak hadronic decay modes as their decay products 
have a much higher total mass. The single particle levels inside a cavity (as for 
the MIT bag model) or for ordinary nuclei or hypernuclei show the same order 
of levels for the lowest eigenstates. First, there is a lsi/ 2 shell, then the IP3/2 
and the lpi/2 shells follow. Due to relativistic effects, the spin-orbit splitting 
is quite sizable for nucleons. As the quarks are much lighter and relativistic 
effects are even more pronounced, the spin-orbit splitting for quarks is at the 
order of 100 MeV for very light bags, i.e. on a similar scale as the splitting 
between the s and p shell. One can put 6 quarks in the s-shell due to the color 
degree of freedom, then 12 quarks in the IP3/2 shell and again 6 quarks in the 
lpi/2 shell. The smallest and most pronounced magic numbers for quarks are 
then 6, 18, and 24 (the next one would be already at 42). 

Studying isospin asymmetric systems reveals another important effect. The 
weak nucleon decay by emitting a proton carries away positive charge. Nev- 
ertheless, the neutron does not carry away negative charge if it is not accom- 
panied by a 7r~. But this decay is suppressed by the mass of the pion and 
the phase space of the three body final state. Therefore, a strangelet stable 
against weak nucleon decay is most likely to be negatively charged. 

Let us look now for strangelets which have closed shells for all three quark 
species with a negative charge and a high strangeness fraction as these are the 
most likely candidates. The first magic strangelet is the quark alpha with 6 
quarks of each quark species at A = 6 which has zero charge [22]. The magic 
strangelets with a high strangeness fraction and a negative charge are then at 
A = 10, Z = —4 (with 6 up, 6 down and 18 strange quarks), A = 12, Z = —6 
(with 6 up, 6 down and 24 strange quarks), A = 14, Z = —8 (with 6 up, 18 
down and 18 strange quarks), and A = 16, Z = —10 (with 6 up, 18 down and 
24 strange quarks). One sees a correlation, that adding two units of baryon 
number decreases the charge by two. These strangelets have a rather high and 
negative charge fraction of Z/A ps —0.5 very similar to an antideuteron but 
with a much higher mass and charge! These strangelets constitute a valley of 
stability which is due to pronounced shell effects. 

This picture holds, i.e. these candidates remain, also within an explicit calcu- 
lation using the MIT bag model with shell mode filling [92] . We calculated the 
masses of strangelets with all possible combinations of up, down and strange 
quarks up to a baryon number of A = 30. Then we look for possible strong 
decays as the emission of baryons (p,n,A, £~, S + , S~, S°, fi~) and mesons (pi- 
ons and kaons) by calculating the mass difference between the strangelet and 
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its possible decay products. For the strong interactions, we also allow for mul- 
tiple hadron emission, like the strong decay of a strangelet via a neutron and 
a pion, and the complete evaporation to hadrons. For example, the strong 
proton decay Q' — > Q + p is checked by 

M(A,S,Z) < M(A- 1,S,Z- 1) +m p (91) 

where M(A, S, Z) stands for the mass of the strangelet for a given baryon 
number, strangeness and charge. Afterwards we check for weak hadronic decay, 
the single emission of baryons and mesons within the same procedure simply 
by changing one unit of strangeness in the final products. The weak proton 
decay Q' — > Q + p is now checked by 

M(A,S,Z) < M(A- 1,S± 1,Z- 1) +m p (92) 

where we allow for both strangeness changing processes of AS = ±1. This 
calculation has been done for several bag parameters. We choose a strange 
quark mass of m s = 150 MeV if not otherwise stated. The value of B 1 ^ = 145 
MeV and m s = 280 MeV is taken from the original MIT bag model fit to the 
hadron masses. 

The candidates which are stable against strong and weak hadronic decay are 
plotted in Fig. 26 in a scatter plot as a function of their baryon number and 
charge fraction. In all the parametrizations shown, we find the candidates at 
A = 10 with Z = -4, at A = 12 with Z = -6, and at A = 16 with Z = -10. 
We do not find any candidates for a bag parameter of B 1 / 4 = 180 MeV or 
higher as strange quark matter starts to get unstable. 

As expected and outlined before, the main long-lived strangelets stable against 
strong and weak hadronic decay are lying in the valley of stability and are 
highly negatively charged. This finding is contrary to the common belief that 
strangelets have a small positive charge and will have serious impact on present 
heavy ion searches for strange matter. These long-lived candidates are still 
unstable against weak semileptonic decay (emission of electrons and antineu- 
trinos) and thus will live on a timescale of maybe 10~ 4 sec [93]. 

One can also ask about the much richer spectrum of short-lived candidates 
which do decay by weak hadronic processes and thus probably do live only 
as short as the hyperons (r pa 10~ 10 s). MEMOs can decay weakly on the 
timescale of the free hyperon weak decay and are thus also belonging to the 
short-lived candidates. MEMOs have quite distinct properties, they can be 
negatively charged while carrying a positive baryon number due to the neg- 
atively charged hyperons, the S~ and the [92]. Light candidates are the 
combinations {2n, 2A, 25"}, {2p, 2A, 2S }, {2A, 22°, 25"}, HO jHe and AA J He 
discussed in section 3. 
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Fig. 26. The charge fraction Z/A for long-lived strangelets, which are stable against 
strong and weak hadronic decay, for different choices of the bag parameter. The 
case for the original MIT bag model parameters (B 1 ^ = 145 MeV, m s = 280 MeV) 
is also plotted. 



MEMOs thus compete with short-lived strangelets as they are of similar 
strangeness content. We calculated light MEMOs up to a closed p-shell and 
checked for metastability (strong decay). We analyzed the strangelet candi- 
dates without the weak hadronic decay, i.e. allowing for the strong decay only. 
The short-lived candidates for MEMOs and strangelets for a bag constant 
of B 1 ^ = 160 MeV are shown in Fig. 27 in a scatter plot as a function of 
strangeness fraction f s , charge fraction Z/A and baryon number A. 

As can be seen, there are many more short-lived candidates than long-lived. 
Light MEMOs can have very unusual charge fractions between Z/A = ±0.6 
indicating a rich structure of strange hadronic matter. Strangelet candidates 
also cover a wide range of charge fraction but are mainly located at negative 
charge. This comes from the strong decay which shifts strangelets to higher 
strangeness fraction and to negative charge. There are MEMOs and strangelets 
with the same strangeness content and baryon number. Here, the energetically 
least favourable object can decay into the other via strong interactions. A 
strangelet created in a quark gluon plasma can then possibly decay into a 
MEMO. Or vice versa, MEMOs can coalesce from the hot and hyperon-rich 
zone of a relativistic heavy ion collision first and then they form a strangelet. 
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Fig. 27. The strangeness per baryon f s (lower part) and the charge fraction Z/A 
(upper part) as a function of the baryon number Ab for short-lived strangelets 
(dots) and unstable strangelets (open circles) for a bag constant of f? 1 / 4 = 160 
MeV. The hadronic counterparts, MEMOs, are shown by crosses. 

Presently, there are only experiments designed to look for long-lived compos- 
ites with a lifetime of r > 50 ns. Designing an experiment for short-lived 
composites is challenging but planned for future colliders and can reveal the 
possibly rich structure of strange matter. 
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5 STRANGE MATTER IN NEUTRON STARS 



Originally (strange) quark matter in bulk was thought to exist only in the 
interior of neutron stars where the pressure is high enough so that the neutron 
matter melts into its quark substructure [18,19,15]. At least in the cores of 
neutron stars (where the density rises up to the order of 10 times normal 
nuclear density) it is not very likely that matter consists of individual hadrons. 
For an overview of the structure of neutron stars see [2,94,95]. 

On the other hand it is also known that the pure 'neutron' matter is not 
really a nuclear matter state made solely out of neutrons, but at least at higher 
densities consists also of a considerable amount of protons as well as hyperons. 
The extrapolation of strange hadronic matter to higher densities in fact has 
influence on the constituents of 'neutron' matter. The gross structure of a 
neutron star like its mass M and radius R is influenced by the composition of its 
stellar material. This holds especially in the case of the existence of strangeness 
bearing "exotic" components like hyperons or strange quark matter which may 
significantly change the characteristic mass-radius (MR) relation of the star. 

In this final section we therefore address briefly on how the two different forms 
of strange matter might influence the properties of neutron stars. This section, 
however, can only give an impression of its implication on neutron star physics. 
It is based primarily on some own work on this interesting issue and can thus 
provide only an idea of possible effects. The interested reader should therefore 
also consult the cited literature for a much wider range of topics. 

5.1 Strange quark matter stars 

In section 2 we have described the equation of state of cold strange quark 
matter by means of a simple noninteracting Fermi gas of up, down and strange 
quarks. In the following we adopt a somewhat modified model which was 
developed quite recently, the effective mass bag model [96]: 

In condensed matter as well as in nuclear physics medium effects play an im- 
portant role. One of the most prominent medium effects are effective masses 
generated by the (average) interaction of the particles with the system. For 
example in the case of a gluon gas at high temperature the consideration of 
an effective mass for the gluon within the ideal gas approximation leads to 
an excellent description of the equation of state found in lattice calculations 
[97,98]. If one considers now analogously the case of a quark gas at zero tem- 
perature the situation is as follows: The quarks are considered now as dressed 
quasi-particles which acquire an effective mass by the interaction with the 
other quarks of the dense system. The effective quark masses are obtained 
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from the zero momentum limit of the quark dispersion relation, which follows 
from the so called hard dense loop (HDL) approximation of the one-loop quark 
self-energy at finite chemical potential [99,96]. They are given by 



for the light quarks (i.e. w,<i-quarks) and 



for massive quarks with current quark mass m s (s-quarks). g denotes the QCD 
coupling constant and is treated as a further phenomenological parameter as 
in the MIT model. As can be readily expected from the formulas (93) and 
(94), the effective masses increase with the coupling constant g and the quark 
chemical potential /i. These effective masses are used in the ideal Fermi gas 
EOS at temperature T = 0. Similar to the expressions stated in section 2, 
the pressure, the particle density p and energy density e of this quasiparticle 
Fermi gas take the following form 
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Here d denotes the degree of degeneracy (e.g. d — 6n/ for rif flavors). Up to the 
additional function B*, which can be regarded as a /z-dependent bag constant, 
these are the ideal Fermi gas formulas at temperature T = for quasiparticles 
of mass m* and chemical potential /i. Due to the /^-dependence of m* (//) the 
introduction of the function £?*(//) is necessary in order to maintain thermo- 
dynamic self-consistency on the one hand and to receive the final expressions 
(as given above) for the particle and energy density on the other hand, having 
the exact form as a noninteracting, but massive Fermi gas [96]. The MIT bag 
constant B still enters as B*(fj, = 0) and is thus a further phenomenologi- 
cal parameter of the model. When inspecting this improved equation of state 
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for strange matter numerically, the overall energy per baryon number E/A 
increases for increasing coupling constant g [96]. This behavior one would in- 
tuitively expect as the masses of the quarks do increase with the coupling 
constant. In return, strange matter thus becomes less energetically favorable 
for realistic QCD coupling constants a s — ^. (One should be aware, however, 
that the value of the 'free' bag parameter is still far from being settled. For 
example, the results of Peshier et al. [98] suggest, that, by including effective 
medium masses in describing the hot gluon plasma and compare it to lattice 
calculations, the employed free bag parameter B^ 4 m 0.7T C might still be con- 
siderably smaller than even the already very low original MIT bag constant 
of Bo 174 =145 MeV. In this sense, all the speculative conclusions of section 2 
are all still valid.) 

We now show how a neutron star would look like under the assumption that 
it entirely consists of cold electrically charge neutral SQM in equilibrium with 
respect to the weak interactions. 

This requires the inclusion of electrons into the model. The thermodynamics 
of electrons at T = can be described by a relativistic Fermi gas 



^ (98) 



re 3vr 2 

£e = 4^ (99) 

* = S- (1 ° 0) 

In return, one has four chemical potentials (fi u , fid, fi s , A*e) which are related 
by the chemical equilibrium between the quark flavors and the leptons. The 
basic weak reactions are given by 



d — >u + e~ + v e - (101) 
s — >u + e~ + v e -. (102) 

The equilibration of flavors is provided by 

s + u — >d + u. (103) 

Hence, the four chemical potentials are reduced to two independent ones 

fi = fi s =fi d and fi u = fi- fi e . (104) 

Finally, the overall condition of electrically charge neutrality ('neutron' stars 
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are locally charge neutral) 

2/3p u - 1/3 (p d + p s ) - p e = 



(105) 



just leaves one independent chemical potential, say p. Therefore, the EOS can 
be written as a function of p only: 



Pb = (p u + Pd + p s )/3, (106) 
e = e u + e d + e s + e e + B , (107) 
P = Pu+Pd+Ps+Pe- B . (108) 

In the following we assume a cold, static, spherical star. It is described by 
the solutions of the famous Oppenheimer-Volkoff-Tolman (OVT) equations of 
hydrostatic equilibrium: 



dp(r) e(r)m(r) ( p{j')\ 47rp(r)r 3 
dr 

dm{r) 
dr 





= 4vrr 2 e(r) . (109) 



They follow from general relativity [100] which one has to apply due to highly 
concentrated matter and therefore curved space-time. The solutions of the 
OVT equations are the pressure p(r) and mass m(r) inside a sphere of radius 
r. The total radius R of the star is determined by the condition p(R) = 
while the total gravitating mass of the star is given by 



R 

M 

r=0 



= J Anr 2 e(r)dr. (110) 



The OVT equation can be solved specifying the central energy density e c = 
e(r = 0) and the EOS in the form p = p(e). 

In Fig. 28 the MR relation of such a hypothetical star is depicted. The nu- 
merical calculation of p = p(e) from (107) and (108) shows that there is no 
noteworthy influence by a change of the coupling constant g. The EOS gets 
slightly softer for high values of the coupling constant (e.g. g = 4) but there is 
no remarkable impact on R and M (see Fig. 28). The maximum mass of such 
a star lies in the range of 1.6-2 times the solar mass and is thus compatible 
with the predictions of normal (nonexotic) neutron stars. One should also note 
one special feature inherent to the description of pure quark matter stars. As 
they are inherent self-bound as a possible consequence of QCD (which just 
corresponds to the speculation about the existence of stable quark matter), 
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Fig. 28. Mass M/M & (in units of solar mass) versus radius R of a SQM star, 
m s = 150MeV, 5 1/4 = U5MeV. 

for smaller masses M ^ M sun , the quark matter stars do become also smaller. 
This is different, however, for normal neutron stars, where their size typically 
increases for lower masses due to the very diffuse (i.e. low density) surface. 

Furthermore, Fig. 29 shows the change of the strangeness fraction versus the 
baryon density for different values of g. Throughout the interior the quark 
matter is indeed strange, with a net strangeness content of f s ~ 0.7 — 0.9. 
Therefore the quark matter star is in fact a strange quark matter star. 

Note, however, that the illustrative results are only valid under the assumption 
of a pure SQM star. Although there is no change in mass and radius of the star, 
there is, nevertheless, an increase of the energy E/A at a given radius inside 
the star due to medium effects. Owing to the increase of E/A in the entire star 
a phase transition to hadronic matter will take place at a smaller radius in 
the interior of the star [94,10f]. The 'neutron' star would then have the form 
of a so called hybrid star, i.e. a star which is made of baryonic matter in the 
outer region, but with a quark matter core in the deep interior. It would be 
interesting to see how the presented improved equation of state will affect the 
picture of such a hybrid star. Such a study has recently been carried out [102]. 
The deconfinement phase transition from hadronic matter to the SQM phase 
is constructed according to a construction proposed by Glendenning [103,95]. 
One here requires the weaker condition of global charge neutrality instead 
of assuming charge neutrality in either phase. The latter assumption would 
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Fig. 29. Strangeness fraction f s versus p/po in a SQM star, m s = 150MeV, 
B 1 ^ = U5MeV and p = 0.17/m- 3 . 

have the drastic consequence of strictly excluding a spatial extended region of 
possible mixed phase inside the star due to a resulting constant pressure in the 
mixed phase. In the mixed phase, however, the fraction of the quark matter 
phase on the one side, and the fraction of the hadronic matter phase on the 
other side, might both be (oppositely) charged if only their volume proportion 
X is always chosen to fulfill the condition of global charge neutrality 

X q QP + (l-x)qHP = 0. (Ill) 

Here qQp and qnp denotes the charge density of each phase. (This argument 
of Glendenning is very similar in its spirit to that one given in section 3 
concerning the idea of strangeness separation during the mixed phase in the 
QGP phase transition.) Since the pressure inside a star must fall monotonically 
from its interior to its surface, a constant pressure mixed phase could not exist 
over a finite region inside the star. However, with the above description (which 
is the thermodynamic correct one) one finds that the pressure varies smoothly 
and continuously with the proportion of phases in equilibrium [103,95] leading 
to a mixed phase of finite radial extent inside the star. 

In Fig. 30 the schematic view of possible hybrid stars of (an assumed and 
fixed) canonical mass of M = 1.4 M is shown for different and increasing 
coupling constant g. The bag constant was chosen so that a pure quark matter 
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Fig. 30. Radii and inner structure of a hybrid star of mass M = 1.4M 
(m s = !50MeV, B^ 4 = 165MeV) [104]. 

star cannot exist even for g = 0. The hadronic phase is described within the 
Walecka model including (some) hyperon degrees of freedom. (Details will not 
be given here and we refer to [102].) One finds that already a small coupling 
constant of g — 1.5 (a s « 0.18) is able to halve the radius of the pure SQM 
core (denoted as QP) of initially R ~ 6 km (with neglected medium effects, 
Fig. 30a) to R « 3 km (Fig. 30c). Already at g = 2 (a s « 0.32) the pure SQM 
core has completely vanished (Fig. 30d), while the pure hadron phase (HP) 
has grown only moderately. SQM is now only present in the mixed phase 
(MP) which dominates the interior structure of the star. For g 3.5 a phase 
transition to SQM and thus the occurrence of the mixed phase is completely 
suppressed. 



5.2 Strange hadronic matter inside neutron stars 

As shown in the previous section, strange quark matter can appear in the core 
of a neutron star. Here we will discuss the scenario that strange hadrons, as 
hyperons and kaons, can appear in the interior of a neutron star. Indeed, it 
was shown by Glendenning that hyperons [105] appear at a moderate density 
of about 2 3 times normal nuclear matter density po — 0.15 fm -3 . These 
new species influence the properties of the equation of state of matter and the 
global properties of neutron stars. There are so many hyperons in the neutron 
star that the whole object is more appropriately dubbed a giant hypernucleus. 
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Hyperons considerably soften the EOS and reduce the maximum mass of a 
neutron star. 



The appearance of another form of hadrons with strangeness, kaon conden- 
sates, was discussed also in many papers [106-110]. Chiral perturbation theory 
(ChPT) gives a rather robust prediction of the onset of antikaon condensation 
at p ~ (3 — 4)po [107] taking into account in-medium modifications of the 
antikaon energy with density. Antikaon condensation will soften the EOS and 
reduce the maximum mass of a neutron star similar to the case of hyperons. 
This will allow for the existence of low-mass black holes and implications have 
been discussed in [111]. In [112,113] it was criticized that effects nonlinear in 
density were not taken into account which will shift the appearance of a kaon 
condensed phase to higher density. Moreover, hyperons were not considered 
in this approach. 

Rather recently, antikaons and hyperons were considered on the same footing, 
where it was found that hyperons shift the onset of antikaon condensation to 
higher density [108] or that it is very unlikely that it appears at all due to the 
rather strong hyperon-hyperon interactions [109]. 

In the following we will discuss the appearance of hyperons and kaons in 
neutron star matter. 

The equation of state for neutron star matter is derived by standard methods 
(see e.g. ref. [105] for the RMF approach without the hidden strange meson 
fields). The equations of motion for the meson fields in uniform matter at rest 
are given by 

d 

da B 

2 * ( B ) 

m at cT =2^g a *BPs 

B 

mlV + dV Q 3 = J29«BP { v ] 

B 

m lRo,o = 9pBTipv ] 

B 

m^o = EWy } • (H2) 



where ps and py denote the scalar and vector densities, respectively. The 
equations can be solved for a given total baryon density ps and charge density 
p c including the contributions from the free electrons and muons 



Pb = J2pv 



(B) 



65 



Pc = Y, VbPv + QiPi = 

B l=e,fi 



(113) 



where qi stands for the electric charge of a species i. In /^-equilibrium the 
chemical potentials of the particles are related to each other by 

Pi = k ■ fx B + ft ■ Pe (114) 

where &j is the baryon number of a species i. This means that all reactions 
which conserve charge and baryon number are allowed, as e.g. 

n + n^A + n , A + A -> S~ + p , ... (115) 

Since we consider neutron stars on a long time scale, the strangeness quantum 
number is not constrained and the net strangeness is determined by the con- 
dition of /^-equilibrium. The above equations fix the fields and the equilibrium 
composition of neutron star matter. 

Fig. 31 shows the composition of neutron star matter for the parameter set 
TM1 with hyperons including the hyperon-hyperon interactions. Up to the 
maximum density considered here all effective masses remain positive and no 
instability occurs. The proton fraction has a plateau at (2 — 4)p and exceeds 
11% which allows for the direct URCA process and a rapid cooling of a neu- 
tron star. Hyperons, first A's and E~'s, appear at 2po, then S~'s are populated 
already at 3po- The number of electrons and muons has a maximum here and 
decreases at higher densities, i.e. the electrochemical potential decreases at 
high densities. The fractions of all baryons show a tendency towards satu- 
ration, they asymptotically reach similar values corresponding to spin-isospin 
and hypercharge-saturated matter. Hence, a neutron star is more likely a giant 
hypernucleus [105]! 

In the following we adopt the meson-exchange picture for the KN-interaction 
simply because we use it also for parametrizing the baryon interactions. We 
start from the following Lagrangian 

L|T = D;KD»K - m\KK - g aK m K KK<j - g a * K m K KKa* (116) 

with the covariant derivative 

£> M = <9 M + iguKVp + igpKrR^ + ig^^ ■ (117) 

The scalar fields essentially decrease the kaon mass, while the vector fields will 
increase (decrease) the energy of the kaon (antikaon) in the dense medium. 
The scalar coupling constants are fixed by the s-wave KN-scattering lengths. 



66 




Fig. 31. The composition of neutron star matter with hyperons which appear abun- 
dantly in the dense interior (from Ref. [109]). 

The coupling constants to the vector mesons are chosen from SU(3)-relations. 
The onset of s-wave kaon condensation is now determined by the condition 

-He = V-K- = WK~(k = 0). 

The density dependence of the K and K effective energies is displayed in Fig. 
32. The energy of the kaon is first increasing in accordance with the low density 
theorem. The energy of the antikaon is decreasing steadily at low densities. 
With the appearance of hyperons the situation changes dramatically. The 
potential induced by the 0-field cancels the contribution coming from the u- 
meson. Hence, at a certain density the energies of the kaons and antikaons 
become equal to the kaon (antikaon) effective mass, i.e. the curves for kaons 
and antikaons are crossing at a sufficiently high density. At higher densities 
the energy of the kaon gets even lower than that of the antikaon! Since the 
electrochemical potential never reaches values above 160 MeV here antikaon 
condensation does not occur at all. The possibility of antikaon condensation 
was checked for several parameter sets and found that at least 100 MeV are 
missing for the onset of kaon condensation [109]. This is in contrast to previous 
calculations disregarding hyperons [107] but in line with the findings in [108] 
where it was seen that hyperons shift the critical density for kaon condensation 
to higher density. 
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Fig. 32. The effective energy of the kaon and the antikaon and the electrochemical 
potential. Kaon condensation does not occur over the whole density region consid- 
ered (from Ref. [109]). 

These illustrations clearly demonstrate the 'richness' of possible scenarios for 
the interior structure of neutron stars. In return it represents a very interesting 
and still open and lively field to pin down the possible structure of neutron 
stars by its observed properties and phenomena. 
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